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1.0 INTRODUCTION 


The early development phase of the Space Shuttle Main Engine (SSME) was 
plagued by a number of problems In which components failed virtually without 
warning. Frequently, these components were in one of the two high pressure tur- 
bopumps which are a part of each engine. Turbopumps, like any rotating machi- 
nery with high speed, light weight rotors have historically been subject tc 
vibration problems that occasionally lead to fires and even violent destruction 
of the hardware. Because of the extreme complexity of turbopump dynamic beha- 
vior, the process of gaining a thorough understanding of the causes has been 
difficult. Testing the elaborate hardware of the SSME has also been difficult 
and extremely expensive. 

The analysis of the dynamics of the SSME and particularly oi the high 
pressure turbopumps presents the analyst with a challenging problem. 
Nevertheless, to gain an understanding of the vibration and whirl problems of 
the SSME, one must develop and exercise analysis and simulation tools. The pri- 
mary tools developed for this analysis consist of a nonlinear, hybrid computer 
simulation and a stability model based on a linearized version of the same 
equations of motion as used in the nonlinear model. 

The purpose of this study effort is to continue work begun under an earlier 
study contract (NAS8-34924) , assisting NASA/MSFC to analyze causes and fixes to 
vibration and subsynchronous whirl problems encountered in the SSME tur- 
bomachinery. Because the nonlinear and linearized models of the turbopumps play 
such an important role in the analysis process, the main emphasis of our work 
has been concentrated on the verification and improvement of these tools. It 
has been the goal of our work to validate the equations of motion used in the 
models including the assumptions upon which they were based and to assure that 
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the linear stability model and the hybrid model predict consistent behavior in 
the regimes where they should. Several tasks were identified at the outset of 
this study: 1. Verification of the SSME rotordynamics simulation and the de- 

velopment of enhancements. 2. Assistance in the review and identification of 
potential mechanisms of subsynchronous whirl and the development and verifica- 
tion of fixes. The primary whirl mechanisms were identified under the earlier 
effort and since an investigation of fixes must await a fully verified simula- 
tion, the primary effort of this study was devoted to the first task. 

The thorough checkout of the simulation has been accomplished and a few 
minor discrepancies have been discovered. These will be explained in detail 
within this report with recommendations made for their correction and an 
assessment of their Impact upon previous results. A filtering technique has 
been developed to enhance the speed ip which results from the simulation can be 
analyzed for frequency content of the signals. Its development and results will 
also be reported here. 
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2.0 SIMULATION VERIFICATION 


During this study our efforts were concentrated on verifying the SSME rotor 
dynamics hybrid simulation and confirming that the linear results agree with 
those of the linear stability model derived from the same set of equations. The 
outputs of the linear stability model are eigenvalues and eigenvectors. The 
eigenvalues reveal the stability properties and the vibration frequencies that 
will be present in the response of the linearized system. The eigenvalues 
should approximate the nonlinear responses as well. The eigenvectors represent 
the normalized amplitudes of motion at the frequency of the corresponding eigen- 
value. Both eigenvalues and eigenvectors are functions of all the model parame- 
ters. The system response is a superposition of all the eigenvectors at 
amplitudes determined by initial conditions and the disturbances acting on the 
system. 

Our verification method consists of making a detailed audit of the simula- 
tion code and comparison with the derived equations of motion and of matching 
frequencies and damping ratios between comparable cases run on the hybrid simu- 
lation and the linear stability model. Agreement between model results over 
wide variations of each of the system parameters gives us a very high confidence 
level (approximately 99 %) that the two models agree both in implementation of 
the model equations and incorporation of an extensive data set. 

Our primary verification technique is to make the hybrid simulation linear 
by eliminating some disturbances and setting all the nonlinearities (bearing 
deadbands, sideloads, squeeze film damping, unbalance forces and internal 
Coulomb friction) to zero. Then we observe the transient response. In this 
configuration, the stability and response characteristics of the hybrid should 
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duplicate the results from the stability model in frequency and decay or growth 
rate of the response amplitudes. 

Since responses of the linearized hybrid simulation contain multiple fre- 
quencies, it is difficult to compare responses to predicted frequencies unless a 
single frequency is dominant. When the system is unstable the instability 
amplitude grows to dominate the other responses at the characteristic frequency. 
The amplitudes approach zero asymptotically in a stable case. The cross cou- 
pling terms of the whirl drivers are used to set up conditions which the stabil- 
ity model predicts to be unstable. The results are examined to assure that the 
hybrid simulation shows the instability at the proper frequency as determined by 
the whirl velocity, the proper direction of whirl, and the proper growth rate. 

The force coefficients which define the force models used in the simulation 
and the stability models vary with the rotor speed (PHIDX). This variation is 
modeled by a polynomial curvefit relating the coefficient to a polynomial func- 
tion of PHIDX. The seals are modeled as quadratic functions of rotor speed and 
in most cases only the term coming from PHIDX squared is nonzero. We studied 
the cross coupling coefficients and their effects on whirl stability and com- 
pared the hybrid arid linear models. 

The third cross coupling coefficients at the Alford force/turbine 
interstage (SQA2), hot gas (SQ12), and preburner pump (SQ22) seals; cross 
coupling at the impel! er/dif fuser interaction point (SQ32), the Alford cross 
coupling multiplicative factor (BETA), and the rotor internal viscous damping 
(TWOZETR) are the stability parameters. For convenience, in this report they 
will be referred to as SQA, SQ1, SQ2, SQ3, BETA, and ZER. Figure 1 shows the 
positions of the whirl drivers on the turbopump model. 
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The methods we use to analyze the hybrid simulation results that are sent 
to the strip chart recorders are derived and Illustrated In Section 2.1. 
Section 2.2 is a discussion of the error due to machine tolerance or operator 
limitations that one can expect in the final results. Investigation of discre- 
pancies during the verification process revealed data input differences and 
simulation coding errors that are described in Sections 2.3 and 2.4. Section 
2.4 also includes an error in the equations of motion that the simulation is 
based on. Our study of the simulation responses indicated other problems that 
are probably due to the hardware peculiarities but were not verified as such. 
These responses o.id implications are discussed in Section 2.5. Examples of hy- 
brid simulation and stability model results are in Appendix B. 
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2.1 DATA ANALYSIS METHODS 


Our analysis of hybrid simulation results concentrates mainly on two para- 
meters: the primary response frequency (Omega) and the damping coefficient of 

the response amplitudes (Alpha). Omega is derived from simulation variables VHP 
and VWT (pump and turbine end whirl orbit angular velocities) which are output 
to strip chart recorder channels 23 and 24. Alpha is calculated from rotor 
displacements in either the Y or the Z direction. We chose displacements at the 
Alford Seal location (DEAY or DEAZ) because these channels (11 and 12) are the 
most reliable pair in terms of pen and recorder behavior. Before any calcula- 
tions are made the direction of whirl is noted. Whirl mechanisms studied in 
this report are forward whirl producers; i.e. the rotor angular velocity (PHIDX) 
and whirl orbit angular velocity are parallel. Any deviation would indicate a 
backward whirl mechanism or a simulation error. The whirl in Figure 2a is for- 
ward, and that of the last part of Figure 2b is backward. 

The magnitude of the whirl orbital velocity was originally approximated 
from the arithmetic mean of the maximum and minimum displacements of VWP or VWT. 
Some simple calculations indicate that the geometric mean produces more accurate 
results. This is demonstrated below: 



FIGURE 3 
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Assuming an elliptical whirl orbit at steady state (Alpha small) as shown in 
Figure 3 and reorienting our coordinate axis to place the ellipse semi-major 
axis along z, we define: 


v w 


y = Ay sin (tot) 
z = A z cos (tot) 
y = Ay to cos (tot) 
z = -A z to sin (tot) 


y,z represent displacement 
to = vibration frequency 

( 1 ) 


Substituting Equations 1 into the hybrid simulation defined whirl velocity, 


yz - yz 


“Ay A z to 

Ay 2 sin 2 (tot) + A z 2 cos 2 (tot) 

The maximum value of VW occurs when y=0, 
“Ay 

vw max = 03 • 

Az 

When z=0, VW attains its minimum value 
“Az 

Vtynin = 03 • 

Ay 

Combining Equations 3 and 4, 


( 2 ) 


13) 


(4) 


“ = Vvw max VWmTrT • ^) 

If the arithmetic mean of VW max and VW mi * n is used to approximate the 
unstable whirl velocity, from Equations 3 and 4 we get 
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t 


1 

1/2 (VW max + VWjntn) = 1/2 (x + -) u> 

x 

- A Z 

where x = . (6) 

A y 

The geometric mean produces a better estimate of the actual whirl orbital velo- 
city. 

Analysis of Omega from the strip chart recorder is illustrated in Figure 4a 
using VWT. The distance indicated by arrow 1 is VWmin; arrow 2 is VW m ax* At 
points after the transient vibrations die out (approximately point 3) measure- 
ments of VW m j n and VW max are made in .01 inch increments. The values are 
substituted into Equation 5 to get the whirl velocity. To minimize reading 
errors, we obtain a mean Omega from between five and ten points. 

The derivation of the relationship between Alpha, the damping coefficient of 
the response amplitudes and those amplitudes is as follows: 

We assume 

y = e- at sin (tut - <p ) (7) 


Where: 

y = displacement in y direction 
o = damping coefficient 
to = whirl velocity 
4> = phase angle 

'i 
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Displacement peaks occur at 


Y = P = e-at when sin (at - <j>) = 1. 

The ratio of displacement of any peak. Pi, and peak Pn, the n th peak from Pi is: 


Pn e 


-o(t + (0 ) 


The damping coefficient. Alpha, is defined in terms of the damping ratio, 5 


1 - 


Substituting Equation 9 into Equation 8, 


= e -2ircn . 


Solving for c. 


1 Pi 

In — 

2irn P n 


And substituting into Equation 9, 


a) In 


1 ■ 

77? - (In Pi/PjJr 


The magnitude, R, and the phase angle, <f>, are calculated from their definitions: 


R = Vo+ to 2 


1 - 5 ' 


12 


yrji 


w£ 


0 ) 

<J> = tan -1 • (14) 

a 

Figure 4b is used to illustrate analysis of the damping coefficient. The 
speed of recorder #1 with channels 11 and 12 has been increased to .5 mm/sec and 
the voltage settings set to 5 volts/line, making the peaks larger and more 
distinct for easier measurements. In this case the peak to peak amplitude 
(indicated by arrows) is measured every tenth peak. Twelve amplitudes are 
measured, the last amplitude (Pio) of one ratio becoming the first amplitude 
(Pi) of the next. Twelve ratios are calculated and substituted into Equation 12 
using n=10 and Omega from Figure 4a. The mean of the twelve values is Alpha. 


2.2 ERROR ANALYSIS 


The hybrid simulation verification process inevitably led us to a need to 
determine "how good is good enough." Several possible sources of error in the 
hybrid simulation results were identified and addressed: data reading and sub- 

sequent compounded error, digital to analog conversion error, and zero drift of 
the recorder. In addition, differences between the stability model and hybrid 
simulation results that may be due to numerical integration and the use of 
single precision arithmetic were investigated. 

The error margin was computed for the case labeled BETB to obtain a rough 
estimate of errors generated when reading strip chart recorders and analyzing 
data. The greatest possible errors for each computed variable (as defined in 
Equations 5, 12, 13, 14) are listed in Table 1, columns 3 and 4. Values in 

r\ 

(table 1 were calculated when we were using the arithmetic mean to estimate 
Omega. Consequently, some values differ from elsewhere in this report 

TABLE 1 CASE "BETB" ERROR (BETA = .5) 


Variable 

Mean 

+/- 

Error (%) 

St. Dev. (%) 

A1 pha 

4.69 

1.33 

28.26 

16.32 

Omega (RPM) 

-20095.83 

251.15 

1.25 

1.13 

Magn (RPM) 

-20095.83 

26.30 

1.25 


Phase (DEG) 

89.87 

.04 

4.24 

— 

The errors in 

columns 3 and 4 

include 

reading error (based 

on j^.005 inches) 


and compounded error in results. The standard deviations of Alpha and Omega for 
BETB are listed in column 5 for comparison with the estimated error. As seen in 
Equation 13, Omega dominates the equation for calculating the magnitude. 


Consequently, the calculated standard deviations of magnitude and phase were 
less than . 1 % In this case. 

It is important to note that a displacement of .01 inches in VWP or VWT 
channels used to calculate Omega, represent) 1492 RPM (156.25 rad/sec) when the 
recorder is set to 5 volts/line. When it is set to 2 volts/line, as in case 
BETB, a .01 inches displacement represents 597 RPM (62.5 rad/sec). The width 
of the pen line is approximately .01 inches as are the lines that represent mid- 
point on the recorder paper. It is difficult for the evaluator to discern much 
less than that. Fortunately the whirl velocity curves are characterized by bar 
type appearances that facilitate easy reading. In this case pen line width does 
not affect the value of the whirl velocity as greatly as the point within the 
zero line from which the measurements are read. To minimize errors it is essen- 
tial to calculate the mean of several readings using the midpoint of the center- 
line. It is with such careful analysis that results of Table 1 were obtained. 

Because Alpha is calculated from a ratio of measurements the problem of 
placement of zero does not dominate the error analysis. However, we no longer 
have -the simple bar like appearance of VWP or VWT. Now the width of the pen 
line is more important. Again a mean value is fundamental to the reliability of 
the values derived from strip chart recorder measurements. 

The digital to analog conversion error was checked by sending a specified 
voltage to each channel and verifying pen movement. Initially a need for reca- 
libration was indicated. All data in this report (except Table 1) were taken 
after recalibration and show only small differences between results before re- 
calibration. Subsequent checks indicated conversion error less than 5% on every 
channel with most channels Indicating negligible error. Because they are used 
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for absolute measurements, the channels of primary importance for verification 
runs are numbers 23 and 24 (VWF.VWT), both of which measured less than \% error. 

No zero drift was observed on either channels 23 or 24 during the error 
analysis. In later runs however, as much as .01 inches has been noticed. It is 
wise to check zero drift periodically on channels carrying absolute data infor- 
mation. 

A simulation was developed on Control Dynamics' PDP 11/23 to isolate the 
effects of the integration scheme and single precision mathematics used in the 
hybrid simulation. Modal damping, stiffness and mode shapes calculated by the 
stability model code are input to the digital simulation. The simulation output 
is compared to stability model results run with the same modal input. Any dis- 
crepancies can be credited to the mathematical methods. 

Figures 5a and 5b are examples of the digital simulation graphics output. 
Figure 5a is the displacement amplitude at the Alford/Turbine Seal. Alpha anal- 
ysis is exactly as explained in Section 2.1 for hybrid simulation output. Figure 
5b is a cross plot of displacements at the turbine end of the rotor. To get the 
frequency from this plot, it is necessary to know the time interval over which 
points are plotted. The number of cycles are counted and the ratio, 
cycles/time, converted to radians/second for direct comparison with stability 
model output. 

A sample of results of the PDP simulation is listed in Table 2 along with 
stability model results. The small discrepancies represent error due to inte- 
gration and single precision calculations. At higher frequency whirls, such as 
those seen here, larger differences are expected (see Reference, page 10). 
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TABLE 2 POP SIMULATION RESULTS; SQ3 = -.0071 


Simulation Stability Model 
Alpha 19.71 18.46 

Omega 4117 4102 

In summary we can say that digital to analog conversion errors, zero drift, 
numerical Integration and single precision arithmetic contribute minimal errors 
in results. The error in our estimate of Omega can be as much as 1000 to 2000 
RPM at 5 volts/line or 500 to 1000 RPM at 2 volts/line. Errors in this range 
can easily arise from reading charts at an angle or from zero drift. However, 
the mean error for a particular value is usually much smaller because errors are 
minimized by careful reading of a large data set and frequent checks of zero 
drift. If results are questionable and greater accuracy is required, the opera- 
tor should evaluate the sources of error to determine an error band for the par- 
ticular run or read data from a digital data scan using the data taking option 
on the hybrid. 
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2.3 INPUT DATA DISCREPANCIES 


In May of 1983 a new set of structural modes and frequencies for the LOX 
pump housing were provided by Rocketdyne and input to the SSME turbopump hybrid 
simulation. Because of the large volume of data required to specify the para- 
meters of the models, it is conceivable that discrepancies between input data of 
hybrid and linear stability models could easily occur. To investigate data 
Input discrepancies, we made an extensive comparison between the two sets. 

Table 3 lists the input data discrepancies found. With the exception of 
the sign of the case mode displacements, the stability model data has been 

changed to match that of the simulation values listed in Table 3. Case mode 
displacements used in the stability model must be the negative of those used in 
the hybrid because of the different way in which they are used by these 

programs. The stability model results are different in the fourth and fifth 
decimal place from output generated before the changes were made. 

The bearing dynamics are modeled as a piecewise linear spring with two 
break points in the stiffness curve. These break points correspond to the inner 
and outer deadbands (see Reference, page E-18). When the deadbands are set to 
zero, the hybrid simulation code nullifies the effect of bearing deadbands. On 
the stability model the parameters affected by the deadbands must be set to zero 
manually when deadbands are eliminated. Consequently, KTDB, the bearing stiff- 
ness between the deadbands on the stability model must be changed from its nomi- 
nal value of 77,000 to zero for cases which will be compared to hybrid runs with 

zero deadbands. The stability model results changed in the second decimal place 
after KTDB was set to zero. All stability model results in this report were run 
with KTDB = 0 unless otherwise noted. 
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Table 3 DATA INPUT DISCREPANCIES 


Variable Hybrid 


Second A1 ford curve 
fit coefficient 

Case mode Y displacements 


Case mode Z displacements 


Rotor mode displacements 
P 
T 
A 

B./S3 

51 

52 

Case frequencies 


Simulation Stability Models 


.7344541 

.5261 

.824220 

-.7791 

-.021071 

.038538 

-.671020 

.51947 

-.030111 

.035695 

-.027291 

.00034801 

.041627 

-.01821 

.6063 

-.33281 

-.038184 

.030607 

.3235 

-.17948 

-.18238 

.11765 

.010491 

-.00918561 

-.89823 

.7987 

-.15703 

.14642 

.58229 

-.46759 

-.30329 

.23216 

.38196 

-.26784 

-.1997 

.12995 

.4776 

-.19665 

-.17592 

.14643 

-.63293 

.37902 


-1.58715 

-1.58458 

• 1.45654 

1.45273 

-.786938 

-.787191 

1.85618 

1.85881 

.08023274 

.08024846 

-4.39593 

-4.398109 

284.062 

284.07 

538.217 

538.27 

700.323 

700.36 

1888.030 

1881.21 

1948.350 

1948.50 

2207.410 

2207.58 

2713.520 

2713.70 

2943.290 

2943.54 

3065.380 

3065.62 

3408.060 

3408.31 
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2.4 SIMULATION ERRORS 


A detailed inspection of the simulation code was made to uncover any 
discrepancies that may have accounted for observed differences in behavior 
between hybrid and stability models. Two very minor errors were discovered in 
the simulation code listed in Appendix A. The stiffness coefficient at the 
Alford seal location (KSA) is not declared real causing it to be truncated to an 
integer value. The declaration should occur at or near line 74. The relative 
displacement from centerline between rotor and case with the outer deadband 
removed (DPP) should be a factor in the multiplication of the deadband stiffness 
coefficients according to Reference, page E-18, Equations 40. This multiplica- 
tion should occur in lines 731 - 746 of the code. The truncation of KSA does not 
affect results significantly. Results in this report are run with deadbands set 
to zero. 


Recently Mr. Steve Ryan of the MSFC Systems Dynamics Laboratory Control 
Systems Division Vehicle Stability Branch discovered an error in the rotor 
equations of motion involving the rotor rotational acceleration, ’<j>x. It is 
necessary to rederive the rotor equations of motion to display this error. The 
derivation will start from the kinetic energy expression for the ith rotor ele- 
ment given in Reference, page E-7. 


Ti = - Mj (Xi 2 + Yi 2 + Zi 2 ) + - Iai^xi 2 + “ hi (“yi* 2 + “zi 2 ) 
2 2 2 


“i " ♦xi + ♦yi sin 4>zi 


4*yi cos ♦xi cos ♦zi + ♦zi sin 4»xi 
♦yi sin <j> x i cos $ Z i + Izi cos * X 1 
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To linearize the resulting system of equations, we reduce the kinetic energy to 
second order equation in the assumed small quantities, <|>yj, <j> zl -, <(>yj, <j> z j. The 
resulting kinetic energy is 


Tj = Mj (Xj + Yj + Zj ) + — I a -| ( <f»xi + +yi 4>zi ^ + — Hi ^yi 2 ‘f’ zi 2 ) 

2 2 2 (16) 


The forces and torques acting on this body can be incorporated into the 
model through the virtual work principle. The work done in a virtual displace- 
ment of the body is given by 

-to lot, O.ppT.r' CACC rr.oocd u-p 

6W = (F x -j5^-j + + F z i ) + T x i (<S4>xi + HiHyi) 

+ Tyi 5(j)yi + T z i 6<j) Z1 ' (17) 


Some discussion of these terms needs to be made here. This expression differs 
from Equation 26, page E14 of Reference, in the additional term which 
multiplies Txi. The extra term is carried along because we allow Txi to be 
large while restricting Tyj and T z j to be small. 

The resulting equations of motion for body i are 



•o 

*i 


>xi" 

Mi 

• • 
Yi 

= 

F yi 






Lztl 


/. zi_ 


T •• y 

H Hi = T xi 


Hi fyi + Hi 4>xi Hi + 
Hi Vzi - Hi Hi Hi 


Hi 4>xi ♦zi = Tyi + T x i ♦zi 

= T zi 


(18) 

(IS) 
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OF POOR QUALITY 

Since I X j ^ - T x -| , we can simplify these equations somewhat by using this 
result in the <f>yj and <j> z j equations. We obtain: 

MiXi ■ F X1 * 

= Fyi (20) 

MiZi = F zl 
*ai ♦xi s T xi 

Iti 5>yi + lai ♦xi ♦zi = tyi ( 21 ^ 

Iti ♦zi “ *ai ♦xi lyi = ^zi 

We now make the substitutions shown below: 

x i = X + jxCj 

Yi = Y + Ij *ij y5j 

z i = Z + (22) 

♦xi ” i>x + Ij Y ijx5j 

♦yi = ♦y + Ij Y-jjySj 
♦zi = ♦z + Ij v ijz5j 

To derive the working equations of motion we make two further assumptions. 
First, the rotor modal functions and ¥jj satisfy axial symmetry and planar 
symmetry properties. Second, the forces and torques consist of the driving 
torque T x acting only at the turbine stages, the internal forces and torques 
which define the rotor mode shapes and T-jj, and the Imbalance forces and 
torques. The resulting equations are: 
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« * ‘ - 
OF POOH *f.'fv..f." 


•X^I.F,,- 
.. 1 

Y = - li (Fyi + F uyi ) 
M 

.. 1 r 

z = — li (Fzi + F U zi ) 
M 


( 23 ) 


♦x = wu 

Vy = -<J»x — lz ~ tlx r o 5y + Ii A i (Fzi + F U zi)]/l2 

12 (24) 

Vz = lx — ly ~ Ux r o lz " Ii*i (Fyi + F uz i)3/l2 

*z 


C y = 


5 7 ~ 


• • • i 


~2?r Ar Ky + <j) X r o 4>y - 4>x r l 5z _ 4>x 2?r Or Cz 

-fi-R £y + £•} <(> 2 -J (Fyi + F U yi) 

-2?r fiR Cz + lx r o lz + lx r l I y + lx 2 ?R n R fy 


-J2 2 r £ z + <f> 2 •,* (F z i + Fu Z i) 


(25) 


The forces F yi and F Z1 * are external forces such as bearing and seal forces. 
Fuy-j and Fu zl - are forces due to rotor imbalances. All forces are calculated as 
derived in Reference. 


In summary, the only changes to the rotor equations shown in Reference, 
page E17, Equation 37, is that all the terms on the right hand side of the 
equality which contained 4> x have been eliminated. The resulting change to simu- 
lation dynamics should be minimal, since typical imbalance forces are of order 
(2500) 2 (lO -4 ) compared to (2500) ( .001) , i.e. 600 to 3. Here we have assumed 
typical rotor speeds of 2500 radians/sec, typical imbalances of 10~ lb sec , 
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typical ramp rates of 2S00 radians/sec, and typical values of Ii on the order 
of unity. Thus in the range of speeds of interest, the erroneous terms are less 
than ten percent of the imbalance forces. As a consequence, we feel that the 
effects of these extra terms on predictions derived from the hybrid simulation 
have been negligible., We do recommend correcting the simulation at the earliest 
opportunity. 

After changing the input data of the stability model and taking those 
errors found in the simulation code into account, we still saw discrepancies 
when comparing the whirl orbital frequencies (Omega) and damping coefficients 
(Alpha) predicted by the linear stability model with those observed from linear 
simulation results. Good agreement of stability boundaries. Omega and Alpha was 
observed when the cross coupling coefficients SQA and SQ1 are varied (Table 4, 
column 1). There was also adequate agreement between the two programs for BETA 
and ZER runs. However, variations in SQ2 and SQ3 were producing discrepancies 
far beyond the limit of acceptable error. The whirl frequencies of the hybrid 
simulation for SQ2 = .005 was 13 % larger than that of the stability model. 
Although the frequency at which SQ3 goes unstable was closer to the predicted 
value, the damping coefficient predicted by the stability model was -32.00 but 
observed to be .87 on hybrid. 

After a thorough investigation of the hybrid simulation and unable to fully 
explain the remaining differences, we ran the stability model with the rotor 
modal mass integral (GAMO), a gyroscopic cross coupling coefficient, set to the 
negative of its nominal value of -.261505. This was done late in the study as a 
final attempt to understand the remaining discrepancies and verify that a pre- 
viously identified correction had been made. The final report of Contract 
Number NAS8-34924 written in October 1982 (Reference, Section 2.1), recommends 
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TABLE 4. STABILITY MODEL - HYBRID SIMULATION RESULTS 


Parameter Value 1 
SQ1 = '.0022 
SQ2 = .005 

SQ3 = .0071 

SQA = .0019 

BETA = .5 
ZER = .15 


STABILITY MODEL 
GAMO = -.2615050 


Alpha 2 

Omega 

5.61 

-2068 

.70 

-4169 

-32.00 

-4179 

9.09 

-2068 

11.03 

-2068 

6.39 

-2070 


HYBRID SIMULATION 
GAMO = NOMINAL 


Alpha 

Omega 

.52 

-2011 

5.60 

-4719 

.87 

-4132 

3.08 

-2083 

4.08 

-2004 

8.49 

-2003 


STABILITY MODEL 
GAMO = .2615050 


Alpha 

Omega 

.35 

-1961 

4.53 

-4650 

.1 

-4160 

2.53 

-1960 

3.89 

-1959 

15.02 

-1958 


1. Parameter value chosen is that at which hybrid simulation goes unstable. 

Nominal Values: SQ1 4.51 E-4 

SQ2 1.576 E-3 
SQ3 1.097 E-3 
SQA 9.453 E-4 
BETA 0.0 
ZER 0.005 


2. Alpha is the real part of the stability model unstable root and the damping 
coefficient from the hybrid. 

3. Omega is the imaginary part of the stability model unstable root and the 
whirl frequency from the hybrid. The values here are in rad/sec. 
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that the sign in the equations using the parameter 6AM0 be changed. Under that 
contract we had determined that the sign of the GAMO term as used in the hybrid 
was opposite the sign used in the stability model. The stability model sign was 
verified to be correct by comparison with the derived equations of motion. It 
was later determined that the version of the simulation being used for verifica- 
tion was not the same as the current version used for production runs in which 
the sign of the GAMO term had been corrected. The correction was verified by 
repeating our check cases with the production version of the simulation and 
noting the response. At that time we concluded that the production version of 
the simulation had the proper sign c. - ' the GAMO term. The latest listing of the 
hybrid simulation code (lines 827 - 844) supports the verification checks run at 
the late stages of the earlier contract. Thus, until all other possibilities 
were eliminated, it did not occur to us to reexamine the sign of this term in 
the current verification effort. We were surprised to see that reversing the 
sign of this quantity as used in the hybrid simulation compared to the stability 
model brought the SQ2 comparison cases to good agreement while maintaining or 
improving agreement in all the other comparison cases. 

Not only does this resolve the SQ2 frequency discrepancies to well within 
accepted error values, but it also improves agreement between damping coef- 
ficients, and whirl frequency values of the stability model and those of the 
hybrid simulation for cross coupling at the other seals. As would be expected, 
the stability model is more stable (with one exception) than the simulation when 
the value of GAMO is used in the linear stability model is set to the negative 
of its proper value to agree with the hybrid model. The numerical integration 
scheme used in the hybrid is expected to be somewhat less stable., but results 
for nominal GAMO were not dependably so. 



Three hybrid runs were made with SQ2 = .005 to examine the simulation 
response to positive and negative GAMO. The whirl velocities are listed in 
Table 5. With the stability model results in Table 4, these runs prove convin- 
cingly that GAMO is implemented with the incorrect sign in the current produc- 
tion version of the hybrid simulation. It is possible this occurred when the 
latest set of structural data was installed. 


TABLE 5. GAMO VERIFICATION RUNS (SQ2 = .005) 

GAMO OMEGA (rad/ sec) 

Nominal -4731 

.261505 -4295 

-.261505 -4759 
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2.5 SIMULATION RESPONSES 


During our verification efforts, we noticed several hybrid simulation pecu- 
liarites. The problems were inconsistent indicating that the cause was not in 
the simulation code. The most likely source is hardware although in these inci- 
dents, there were no run time error messages or malfunctions noticed. 
Diagnostics did not indicate any other problems. This section, therefore is 
more for completeness of our report than to point out any simulation errors. 

In August there arose a question concerning the sign of PHIDX, the rotor 
spin angular velocity. The seal coefficients are express'd as polynomial func- 
tions of the shaft spin speed called PHIDX. Since PHIDX is negative for the LOX 
pump the absolute value of PHIDX is used in these polynomials. The cross cou- 
pling coefficients should be multiplied by the negative sign of PHIDX as their 
affect opposes the rotor velocity (see Reference, E-17, Equation 37). We could 
not find the multiplication in our copy of the code. 

As a result of the forms of the equations and the way in which the model 
coefficients are calculated from PHIDX, the hybrid simulation should behave 
identically except for changes in the direction of induced whirls when rotor 
spin direction is reversed. Test cases run in September with all destabilizers 
except SQ3 set to nominal failed to produce the symmetry in whirl velocity 
magnitudes for +/- PHIDX. Testing the other parameters with +/- PHIDX produced 
symmetrical whirl velocities. Relative to results we were observing for SQ2 
stability-hybrid comparison runs, the output was close to the outside of our 
error margin so we turned our efforts to the more alarming discrepancies in SQ2 
that we later found to be due to the 6AM0 error. 


In late December we returned to the PHIDX sign problem. To insure proper 
treatment of the cross coupling coefficients, we checked the response of the 
simulation to reversed shaft speed (PHIDX) and reversed cross coupling coef- 
ficients at each of the seals. All destabilizers (SQ1,SQ2,SQ3,SQA,BETA, and 
ZER) except the parameter of interest were set to zero in order to isolate its 
effect. Four runs at an unstable coefficient value were made for each parame- 
ter; +/- cross coupling coefficient and +/- PHIDX. We were looking for proper 
response in direction and magnitude of the whirl velocity. The whirl velocities 
responded correctly in direction for sign reversals at unstable values of 
SQ1, SQ2, AND SQA. For an unstable value of SQ3 and PHIDX greater than zero the 
whirl velocity did not reverse direction until it had built to amplitudes too 
large to measure. 

Because of the discrepancies in these results, the SQ3 test (with all other 
destabilizers set to zero) was run once more in February (Figure 6). This time 
the whirl velocity responded correctly in direction. Magnitude differences were 
no greater than 1480 RPM. This run is the most reliable. Although our test 
cases did not produce questionable results at the beginning of the September and 
December sessions in which the SQ3 runs were made, other runs on those days 
indicated possible hardware or system software malfunctions. These problems 
that were seemingly circumvented are enough to seed doubt in those results. The 
"cleanness" of the February result provides us with additional reason to disre- 
gard the earlier results. 

Our data indicates that the sign change in PHIDX is being treated correctly 
in the calculation of the simulation parameters from the curve fitting polyno- 
mials. We cannot verify this by inspection of the simulation listing but note 
that it is apparently being done and are content with that. 
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<— Figure 6d 
VWP « -4105 
VWT «= -4254 
SQ3 « -.0085 
PH I OX * 3194.5 
2 v/1 Ine 


Figure 6c — > 

VWP = 4259 
VWT = 38.05 
SQ3 = -.0085 
PHIDX = -3194.5 
5 v/1 ine 


< — Figure 6b 


VWP = 4429 
VWT * 3906 
SQ3 = .0085 
PHIDX = 3194 
2 v/1 i ne 
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Figure 6a 
VWP = -4285 
VWT = -4433 
SQ3 = .0085 
PHIDX = -3194.5 
2 v/1 ine 


VWP VWT 



VWP VWT 


FIGURE 6 SQ3 VERIFICATION RUNS; 
FEBRUARY 1984 


VWP, VWT, and PHIDX expressed in 
rad/sec 
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We Investigated hybrid simulation repeatability In October. This was done 
by rerunning several cases using identical parameters. Excellent results were 
obtained for BETA and SQ3 comparisons." SQ2 results were not as good (Table 
6a). 


To find the stability boundary often requires several run time inputs; the 
simulation is suspended while the operator enters a new parameter value. The 
parameter is changed from unstable to stable values, then to an unstable value 
closer to the boundary. The process is repeated iteratively until the stability 
boundary is found without letting the system go unretrievably out of control. 
What we will call "cold start" runs are those data runs not preceeded by any run 
time inputs. 

For cold starts the SQ2 results are within acceptable error margins (Table 
6b). The data in itself does not indicate whether the problem is due to SQ2 or 
hysteresis. Rerunning the simulation with a positive GAMO may produce more 
conclusive results. From Table 4 it is evident that GAMO does not affect all 
parameters to the same degree, alghough one would expect it would affect all 
runs of one parameter similarly. From these results, we must assume that there 
is no connection between earlier SQ2 hybrid vs. stability model discrepancies 
and the relatively large standard deviation between hybrid runs. 


TABLE 6a COMPARISON OF WHIRL VELOCITY 


Parameter 

Run Date 

OMEGA 

(RPM) 

Run Date 

OMEGA 
• (RPM) 

*D1f 

Beta=3.0 


8/30 

-16644 

10/19 

-16664 

.12 

SQ3=.0071 


9/6 

-39458 

10/20 

-39515 

.15 

SQ2= t 005 


9/6 

-46028 

10/20 

-43354 

6.16 

TABLE 6b 

COMPARISON 

OF VW (RPM) 

; SQ2 = 

.005, COLD STARTS 


11/10 

11/28 

11/28 

Mean 

St. Dev. 

% 


45636 

44853 

45436 

45311 

407 

.90 
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3.0 FILTER 


We have developed a real time adjustable filter that is capable of 
filtering Y or Z displacements at any joint to find the dominant whirl fre- 
quency. The damping coefficient can also be derived from the filter output. 
The filter uses a running Fourier coefficient algorithm. The derivation of the 
algorithm is included in Sec. 3.1. Section 3.2 explains the use of the program 
and presents some examples. The FORTRAN Code is listed in Appendix C. 
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3.1 FOURIER FILTER ALGORITHM 


A sampled function f(t) that is periodic with a period of x = NAt can be 
represented as a Fourier series 


00 2nkt 2irkt 

f(t) = l (a|< cos + b|< sin ) 

k=0 x x 


( 26 ) 


where 


a k = 


N 

I 

n=0 


2irkn 

f (t n ) cos ; t n = n At 

N 


N 

h = l 

n=0 


f(t n ) sin 


2ukn 


For convenience, let 
fn s f(t n ) 

2Tikn 

Ck,n = cos anc * $k,n = sin 


2-irkn 


(27) 


(28) 


At a time t = t + At, 

a k' = n+1 c k,n 

n=0 


b k' = fn+1 s k,n 
n=0 


(29) 


Manipulating Equations 29: 
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I 


N 

a k'Ck,l - b k ’S kj i = I f n +l (C kj n C k ,i - S k>n S k>i ) 

n=0 

N 

= l f n+l c k,n+l 
n=0 

= a k - f 0 Ck,o + f n+l c k,n+l 
N 

a k ' s k,l + b k 'C k> i = I fn+l (C k>n S kj i + S kjll C kj i) ' (30) 


N 

= 1 fn+l S k n+1 

— — rv 


= b k - fo s ko + fn+l s k,n+l 


From Equation 28, 

C k ,0 = 1 C k>n+1 - Ci 

S k ,0 = 0 S kjn+ i = Si 


Substituting Equations 31 into Equations 30, 

a k ’C k ,l - b k ’S kj i = a k - f 0 + f n +l C k> i ^ 

a k 'S k ,l + b k 'C kj i = b k + f n +i S kj i 

Solving for a k 'and b k ' simultaneously, 

a k' = C k) i (a k - f 0 + f n +l C k i) + Sk,i (b k i + f n +l S k i) 

(33) 

b k '= S K> i (a k - f 0 + f n +l C kj i) + C kj i (b k> i + f n +i S kj i) 

In Equation 33, the coefficients a k 'and b k 'are expressed recursively in terms of 
the first and next values of f(t) in the Fourier interval, the immediately pre- 

ceeding coefficients, and the cosine and sine of that frequency for n=l. By 

graphing the coefficients or the magnitude of the two for a k value, one can 

observe the trend for that frequency. 
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3.2 FOURIER FILTER OPERATION 


The Fortran 77 listing in Appendix C is implemented on a Hewlett Packard 
9000. Although we run the filter in series with our digital single precision 
simulation, it is designed to accept data in parallel as a subroutine might. 
There is a descriptive list of variables at the beginning of the computer 
program. 

From Equation 26 we define the frequency interval as the fundamental fre- 
quency'll) of the input function, f{t): 

2it 

u>o = — s Atu (34) 

T 

For a sampling of N simulation points each separated by the simulation time step 
DT, 

x = N • DT 

2tt (35) 

Aw = 

N • DT 

The user selects a frequency interval (FREINT) and inputs the simulation DT. N 
is calculated from Equation 35 to insure that the user's choice of filter fre- 
quency (FILFRE) will be an integral multiple of Aw. Because N is an integer 
and it is essential that fo is exactly N + 1 points from fn+l. Aw is recalcu- 
lated within the program using Equation 35. 

The wave number, k' that will give FILFRE (w) is 

w 

k' = (36) 

Aw 
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Up to 4 other frequencies are chosen on each side of k' to give a frequency 
band. 


The program uses either simulation time, number of sampled points, or 
number of Fourier intervals to signal termination. The program will stop when 
it reaches the first of these conditions. In our program we control resolution 
of graphics output by changing NPTPLT, the number of points per plot. 

Output variable, AK, BK and RMAG are sent to a file in a format acceptable 
to our’ plot programs. We have found that RMAG plots are the most easily ana- 
lyzed for the high freqencies we are interested in. RMAG for several frequen- 
cies can be sent to the analog parts of the hybrid for output to strip chart 
recorders. The result is a real time frequency filter that can indicate domi- 
nant frequencies while the program is running without having to wait for tran- 
sient simulation variables to settle out. The filter frequency can be adjusted 
during a run-time input. 

The program filters nine frequency values; the highest coefficient magnitu- 
des on the plots belong to those frequencies closest to the unstable frequency. 
In this example the simulation on the PDP was run with SQ2 = .005. Figures 7 
through 10 graphically represent filter output at the Preburner Pump End 
Bearings. When the filter frequency is set at 3400 rad/sec with a frequency 
interval of 100 rad/sec, coefficients of frequencies from 3000 to 3800 rad/sec 
are plotted. The largest coefficients occured at 3800 rad/sec with little 
slope, indicating that the unstable whirl frequency is outside this range 
(Figure 7). 


When the filter frequency is set to 4300 rad/sec the coefficients at 4200 
rad/sec (RMAG(4) have the greatest magnitude (Figure 8). Successive runs with 
frequency intervals set to 50 and 25 rad/sec (Figures 9 and 10) limits the 
unstable frequency possibilities to between 4162 and 4188 rad/sec. The effects 
of displacement along the rotor axis can be seen by comparing Figures 10, 11, 
and 12. It is possible to derive the relative mode shapes in this manner. 

The frequency interval, Atu, and number of points required for a Fourier 
interval are inversely related; a decrease in Aco requires a larger number of 
points, N, and a greater sampling time, t, before acceptable data is available. 
The time step between sampled points can be increased, to reduce the number of 
points sampled and therefore the number of data points kept in an array. 
However, a greater sampling time must, still be used. 

For example, in Figure 9 the number of points in the Fourier interval is 
4117, and the frequency interval is 50 rad/sec. The data is good after .127 
seconds. When the frequency interval is reduced to 25 rad/sec, the Fourier 
interval is maintained at 4117 by taking every other displacement point from the 

_5 

simulation, i.e. DT ~ 6.104x10 . It is clear from Figure 10 that the data is 
acceptable after .25 seconds. Although we have actually increased the amount of 
simulation time necessary to filter frequency to a resolution of 25 rad/sec, the 
number of data points in our array is maintained to prevent twice the number of 
calculations. We have taken up to every tenth simulation point with no obser- 
vable differences in graphs. This is not necessarily a limit; but because it 
would require a very long simulation run, we stopped there. 
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The analysis of filter output is illustrated in Figures 13 and 14. These 
plots are output from a test case where a simple sinusoidal function, F(t) = exp 
( .005t)sin( .63t) + exp(-.002t)sin( ,88t) was input to the filter. The peaks are 
measured on the AK coefficient graphs and alpha computed using the same methods 
as used for hybrid simulation displacements. The dominant growth rate with the 
filter set at at .628 rad/sec is indeed .005. The growth rate can also be accu- 
rately derived from the magnitude plots by using the relationship derived here. 


Magg = Ke at o at t 0 (first point in sampling period) (37) 

Mag = Ke + At)] (38) 

Mag 

= e <*At (39) 

Mago 


Mag 

In (40) 

Mago 

o = 

At 


According to Figure 13 the output is not good until after 100 seconds. 

Picking a peak after 100 seconds on Figure 14a we measure its peak to peak 

amplitude (1) and that of another peak (2). By choosing n=l in Equation 12, 

Sec. 2.1, one can get many readings to eliminate possible error due to time step 

choice or poor plotting resolution. In this case, 

412 

In 

.628 393 

a = (41) 

V* 2 - (In — f 
v 393 

(Lengths are measured in .01 inch increments). .63 is the filter frequency which 
can be checked by reading number of cycles per unit time off graph and convert- 
ing to radians. 
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F(T) = e‘ 005T sin(.63T) + e' ,002T sin(.88T) 
NN (no. of points in sampling period) = 100 
DT =1.0 sec 
Filter Frequency = .63 
Frequency Interval = 6.28 x 10"2 

For k = 10, 0MEGAK = .628 
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FIGURE 14 FOURIER FILTER OUTPUT 
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Using Equation 29 and Figure 14b we get the same results by measuring from 
0 to same point in cycle for several places along curve. Using points 1 and 2, 
252 

In 

153.5 

o = = .005 . (42) 

100 
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4.0 SUMMARY 


Under this study we discovered some minor errors, none of which are 
expected to cause significant changes in predicted bearing loads or responses, 
although they do cause the stability boundary to shift upward. For this reason, 
the sign of the GAMO terms should be corrected on the production version as soon 
as possible. Although not as significant as the apparent sign error of GAMO, 
the other errors discovered need to be changed to assure the accuracy of the 
hybrid model and its agreement with the stability model. These include dropping 
those terms in the equations of motion multiplied by the rotor rotational acce- 
leration, declaring the Alford seal stiffness coefficient (KSA) real, and 
including the relative displacement DPP when multiplying the deadband stiffness 
coefficients. 

We found that for reading whirl frequencies from the strip chart recorder, 
an error of 1000 to 2000 RPM, at a strip chart recorder setting of 5 volts/line 
is possible. However, using lower volts/line setting to enlarge recorder out- 
put, meticulous reading of charts, and averaging several readings can be used to 
minimize errors. With these techniques, errors of 240 to 400 RPM when recorder 
is set to 2 volts/line are achievable. 

There are some inconsistent hybrid simulation results that we were unable 
to irrefutably attribute to any one source. The primary example being the SQ3 
cases presented in Section 2.5. We were unable to determine why recordings of ± 
PHIDX In December were unsymmetrical , while those run in February react as 
expected. There have been other instances where Identical input did not produce 
identical output. Our attempts to characterize these inconsistencies did not 
demonstrate anything conclusive, although there are indications that hardware 
malfunctions and/or hysteresis are the causes. 
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The Fourier Filter algorithm we have developed would be a valuable enhan- 
cement to the simulation, allowing for real time frequency analysis. The filter 
uses a recursive Fourier series algorithm to analyze displacement amplitudes at 
different frequencies. Several frequencies could be evaluated with the results 
for each going to a different strip chart recorder channel for analysis. 
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CONCLUSION 


The primary concern of this study has been to define and quantify any prob- 
lems in the SSME hybrid turbopump simulation model. This has been done by a re- 
view of the derivation of the system equations of motion and by a painstaking 
comparison between a linearized stability model and the nonlinear hybrid model 
where both have been set up with the same data and conditions. We have made our 
recommendations of corrections for the four errors we have discovered and we 
feel these corrections will remove the remaining differences between the hybrid 
and the stability models. The problems we have discovered do not invalidate the 
results run in the past but the corrections we have recommended should result in 
closer agreement between the linear and hybrid, nonlinear models. From our com- 
parisons between simulation results, we can say that, with all recommended 
corrections made, the SSME hybrid turbopump simulation correctly models the 
equations of motions, including the known whirl drivers and significant nonli- 
nearities with sufficient accuracy. 
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OUTPUT U)>Kt|<KSt«^l 

OUTPUT |eiEP8C»eKPY*BKP|,e K TY»KKTj 

OuTPuT *PPTO#AKFYl*AEPY|,AKPYJCiNPY^ ' 

0U!PuI*J?m/i£fc2 

OUTPUT AKTYQ#aRty}»* k TY*»aRTyJIAktY< • 




. OUTPUT PHIpPiPHjRT^PHIRA^pMToetPHlPSl^MlPSa 
OUTPUT PHJRJ3 

— iaah^ai:WTa:ai3sa» Uic,..h.,et.,>,i 
SSifcl ES!JSli:?!:li*i !" ,MC, * TtPH ’ e * 11 - ■ - 

OUTPUT {Cl/A;j»AC3#*CY,AC8#ACtEACT 

OUTPUT 6H»*N#0HPRsi/0NrC#8V'NCsc 

0WT?Ul..9PiBIffl?Ei?pT 

OUTPUT JCAV#RP*RT'UPP#UPT#ElP«TUT#CHAXP,EMA*t 
OUTPUT CtP»ClTi(2^CET f C3p«t?T ... 

OUTPUT SKAC*SKAi ( SKAt<SCAC«8CA2 # S0Al i 9CA0«SCAt»SCAi 


(3/ 



V 


XI * 


- .372 •_ 

37*. 
375. 
37*. 
3)7. 
37* . 

- 37$.- 
3*C. 
3*1. 
3*« • 
313. 
3**. 
-215 — 
3*6. 
3*7. 
38*. 
3*S. 
39C. 
-391... 
39c. 
392. 
39* • 
39*. 

396. 

397. 

39* . 

» 395. 

0 *CC • 
*C1. 
*Cc • 
'C3» - 
'C'. 
'C5. 
'Ct. 
'C7. 
Ul. 
'cs._ 

'1C. 

'll. 

' 18 . 

'13. 

15 

'It. 

'17. 

'it. 

'is. 

'5C • 
-J5l._ 

' c c . 

'83. 


OUTPUT 8MC>SKli,SK|g,SQlC'S0l5«8Qlt,tci0'! 

_0LIFut_SK2C<SKt {. S kCSj S02C<SCJS.Ja?f*Jcic<! 
OUTPUT 6»»3c#SK3i,£K3I.S03C**335«Sa32,8C3C<! 

OUTPUT FSlMX 

OtTPuT CrY#CP|#CT7,CT2 

OUTPUT CPPt«CPPl.CPTY,CPTi 

OUTPUT AK6C*AKei,AK8I#ACBC«*CB»,ACB5 

.6LIEuT-£POeY,CPc|2,C108J, CTDB2 

OUTPUT UPPY»UrPl,UhTY#ur*TZ 

OUTPUT .UPIY»UmBi,UKAY«UPA| 

OUTPUT UP|Y»UPl|,Lh2Y,uP5l 

OUTPUT UP3Y*UM3: _ 

OUTPUT OSPEEQ 

-OUTPuI— IX02ETB, JW02E1C 

OUTPUT *PUMX|Oc 

OUTPUT I!PE#8PEj6 - - 

OUTPUT T|KtOk#TcuPP»NPAA#lMA*«SPA*iO 

OUTPUT 01, TF ... 

OUTPUT 0T8PlB,>sC0w 

. flUIEuI_JD0,10,Jj iCD#20i 2 ! 

OUTPUT PH|COX#Ph|Cx»PHIX 

OUTPUT PH|OCYiPklCY#PH|Y *. 

OUTPUT FH|002,Pu|C2,PMI2 

OUTPUT M0CY,x!CfAX|Y,XtD02»Xtl2,x!Z 

OUTPUT ETAC0,CTA0*5TA 

_OUTPuI_£P.l« 8C Alt.* > SCALE, ttH£!M 

OUTPUT TPHIOa 

OUTPUT lOCT8T,OcBT8T/T!HT§t..„ 

OUTPUT I S |CEL 

OUTPUT TH6U0S , _ 

OUTPUT TH$Ul6uk ( $UY,5LI 

-OUTEuI— iH^UAi iLir/SLATj Sui2 


OUTPUT Tp$U0l,SujlK*5U$lY,SlSlV 

OUTPUT TKSLOZieujJK/OUSiYiSLSJV 

OUTPUT TH$tn»Su«JMSu$TY»*C$lV 


OUTPUT ASPC0Yi*KFCB2/AKTD6TiAKfDB2 

OUTPUT ABF 

OUTPUT -tCAPPJ,EtiPPl - 

OUTPUT FCAPTY,Fc*PT2 

OUTPUT KPUPP 



OUTPUT * 

C 

*P.I1eU?I»?5ZJ 

AB I Tc C 151,787 1 

%P£PoO?C_ 

CAUL R3L nREhODjSlRPOD) 

IFURErOCIuPErOoil 

|F(jrE»*OC»EC.I»ob TO 20CI 

CALL_B5L UTOEiSutOEJ 

|F«lT0F)WRlTE«lC8*757) 


c 

C 



r t 


1 


J26. 

*2 7. 
* 25 . 

*2S • 
*3C. 
*31. 

.*32. 

*33. 

*3*. 

*35. 

*36. 

*3* • 

*25. 

*3$. 

**C. 

** 1 . 

**c • 
**3. 

.***. 

**5. 

**(. 

**7. 

**5 * 
**S. 

*5C. 

*51. 

*52. 

*53. 

*5*. 

*55. 

.*56. 

*57. 

*55. 

*5S • 
*6C. 
*61. 

*6 c • . 
*63. 
*6*. 
*65. 
*66. 
*67. 

*66 • 

*6S • 
*7C . 
*71. 

* 72 • 
*73. 

*71 . 

*75 . 
*76. 
*77. 


_c c ;: <mm 

^ 1 CONTINUE “ 

C THIS SECTION pRevlOES THC.ABJl8TT T6_0aVHER DitA DuPlNO.tXteC* I0N_ 
l^iJOAlHggl_Oa_tfl_l*Q l _ 

CALL R$L( lOATHCp.SLOATHEfl) 

m*K0l* IQMHEm 00 TO i*og _ 

co uci 1 1 1 # 36 ~ ~ 

QMrAMn«cAni) . t • 

OA|P|M|)ftDA|(it “ 

H-2-COMlNUC 

CALL wDAC(C«35M£PO>Srn 

1 * 1 COM IN VI _ • ‘ 

KGFAsS^OPASSAl 

IMXC?AS8i0Tt0fxfiPASS)0« TO i*U 
00 U03 | il# j5 ~ 

|% 3 C«KT|NU| _ * 

1% C C 8 M 1 MJE 


-wHwaiWiWWhi 

GO Tg l! 


Cl*UKSO# MUm»UN » oo TO 10 


OAT* WA« OATHEBJC THH T|Ht 


JliC-EgEgin?^UhiJDLi*jfiH» Jurf indcd .i 

lA'l FORPaTUSH OOt CAIA* RONifO • Eaa# JH» #lJ»l?HJUST iNDfOi) *- 
MI2 FORHATjjKlAllH OAjA.fOR BfM|0« .* l?» PXiAHCHAM3X#_ ... 

• WM DESCR|Ptofi >IX«12H H#X« vltUI «)X#ltH Pin* VALUE# *7I#F» «*# 


»Aii*F0fil'AllMj«#fi.«x;3AE;3x#Ei8^»’S#tIf.*» 

__..njA.coMifcy! , . . : .. 

X R IT E ( IT VP 1 1 1 U l ) ffoTfoTNUpffUS 
El^fTl^FjKOFASS^Oj 

*>R ITe * tC.A J.l* IficMO* nup«UN#eUapt}pe 

co.lvo? M»*3I 

J*|C|.NNL|K) 

1**5 COKT 1 Nvi 

9UTFUT iesEPBc#eKFY*exPbeSTy#KXTl. r _ . 

OLYPlT **PYC#aKFyJ#AKPYE»aKPyJUkPY* 

OUlPtT AlfpjCfAKpjbAKP^AKPmAKPE* . _ 

OUlPuT APYY6#ARTYi#AXTV? # AKTY3tAKTY* 

6WlELil-AKI26i0^.tZl«AinZ(i A^l2)|AKIZA 

OLTPlT *'PV#kPZ»xtY#kTZ#xprY#kp9z,kPTv7kptI 
OUYPul Pt'fPzOHf’RS?--. 

OLTPlY OP,flt»GPp,CpT 


<w 



rwrvii 


Ml. 

05 ^ 

Ml. 

Ml. 

Me. 

M3. 

MM 

MS. 

Mt. 

M7. 

M!. 

MS. 

MC. 

Ml. — 

M c. 

M3. 

MM 

MS. 

Me. 

-M7. 

Ml. 

MS. 

5CC. 

SCI. 

SCc. 

SC3. 

. -0M 
* -CS. .. 

> see. 

SC7. 

SCI. 

5C9. 

5|C. 

Sit. 

Sic. 

513. 

SM. 

SIS. 

Sit. 

Si?. 

S|t. 

sn. 

s:c. 

Sc 1 • 

Sic . 
S13. 

sit. 
Sis. 
Sit. 
517. — 

Sic. 

SIS. 

S3C. 


OUTPUT lCAy«RP*RT*0 , 'P»^T*StP#Ttt#tHAlip,tHA)(T 

OLlEut— ClP.CllAClP. Ctl. C3 PaC31 : 

OCTFuT F$ 1 HAX 

OCT Put CPtiCPI/QTTcCTI ._ . __ 

OCTFuT CPPYaCPPj.CpT Y.CPTZ 

OCTPfT AKBC«AKBi.AkB2. ACB0*AC81*AC62 ... . 

OCTFuT CPoeY/CPOB?/CTOBY/CT08T 

BCISuT-CPPYcCKEl.LKTI.UtTl— 

OCTPuT CKOY»CKBi,ChAY«UHAl 

OCTPuT CKlYcUHli.CMlY.Ctll 

OCTFUT CKJY#CMJ2 

OCTPcT IS JCEC.Tu jLe5 

octpuT thjc»sc*»sCyc»scz6 

OCTEul— ItSLA. SLA^i SLAYO* SCAZO , ! 

OCTFuT TKSLSI.SljUaSLSIYC'SLJSZO 

OCTPuI TK$L82>St,ftKiSCStYC*6lS3!Q . ‘ 

OCTPuT THfiLS3.Sus?KASCS3YC<SC6AlO 

lOATuCR * .FAtSc. . I.. 

KOFAsS*C 

THJ AREA BESEPycO FBR RnNItUtlilNO AfTtR A RUNTIME INPUT 
1C CONTINUE * " 

input i if ice I -s- 

NCKRuN*M.'HRUNct 

WjDt2HIS*X8L0H 

0181.0I/1. 

TNOOT'OTOT - - - . 

OY10t*OToa«OTO* 

C1F*1 . •CKPcELP* i(lP/(0PP»(lPll««5 _ . 

C1T*2 »*CKT*ECT* (RT/(0PT«0T)I <, A 

cae.icLt/e*/RE ; 

C3T.ELT/1./RT 

ClPrl.tClPot ... 

ClTci.^ClT*** 

ELOOiP*C3Pc»2 . 


tL00iT*C3P.*l 

CLEJOet.*(Jl! 

CLTtoPT»OT 

. OhKRsO^OKHRtOWtB. ...... ; , 

IFt lS|Ott.fcC«.t » 8 l** 0 I ICIlKlO | Zcl2K*0 I sLSSK.C | IUaMO 

CO lol PYTaKOOC, _. ... ... 

>1 OKKC$C«KMOHtC(K)*OKHClR) 

ll_CONTlNUE 

CALL P5CMEX1T/1J ..." . 

IFllEXlTl CALL IxlT 


. JEIKTAPE.NEiOI- 00-JS.71C 

CALC R81UT*PEM$CTAP) 

IF ( ITAPE )G 0 . TO.?cS . .. . 

|F IT «LT • T IKEON Ifla TO 71C 



v *7 ; & k * 




T INEflNHCUNP ( | 

? 9 JPASStP * 

NPASS'O 

mahfOt ! 

K!APe»«1 

KCWt'.ec^DUKPil , 

KTEKp.KCCPP*| iOt>P 

CAU WCUKTENP* 

c 

E VRItE VEACER9 M tH|8 t !MC-SOlCf 

C 

fcPHEtRIAEEtZBCj^HAtftieAliNHAlUjDlAAuNlP 

*>P|U(» , TAPEa7B1|(A8CT(1)#!*T»11) 

VP|T£jUPp#75E)NCHANtlrAA|Nf’{*#i/OTEf'ii l RUN|0 , 

hR|TE|E' , P*73J)(»jeuTU)#AeU!U l )/|«Uio» 

7 c F 8 prAn:</ 2 iA, 8 i 2 c, 8 /A») 

7 1 FORPaT ( lCAA ) 

Z^J-ee^I^/^ajAVbAN^ELi/fAhOtEj KJiti.SAMfLES/RfCOBQL 

• ?ApplE 6 /T|pE. 1UCE. TlPn^ElA,7j I (WITH OT • • «ElA »7i 

• PVN»IO« * # A % ) 

7 3 feSPATM cHANNiLti'^iorr'TrnfiA'nri^n — 

7 C CONTINUE 

c 

£AUi_E>£klAS 


PHIOjiSO* PH|05I*PNl0A 
i PHjOAA*AesiPHlDij 


AKf Y# CAKgYC*PH|0)fA?fAKPYl#AK5y|«p H |pyA^A«PY3.p 
AKFf.lAKPIC*PHI»;AA*«AKP7l4ARp;J»PHl0xA*AKP73.p 
AKTT*(AMtC 4 PhIi 3 *AsIAMTY 1 *ARtYI*PhIDxA»AKTT 3.P 


JtllllAiLllCiEM 

AKTI«|AKTZC«PH 


pwiOXS9M!AkPYj 

PHIDX$C))"AKPZ| 


IClEHlQX^iAKIll±ARlYl#EHiDX4iAKliaiEmD»SJlIlAmi 

ZC«PHlOxA«IAKTI14lKTriipHl0xA4AKTZJ.pHIOXSCn«AKTZ4 


C KIIMf** lECTiflKJF $TJ!FN|68 CeNftAli'TSFOP.BEARiNOS.ANO THE t W 

c respective carRjeRs ASe to be RBnsIDereo InoepenoenTly; 

~c AK (P.T ) ( Y 47 ) (0*4 ) APE THE CURVT*FJt C8|FF fC TenT 3 Fe« BEaRTnUI# 

C BKJP.THY.ZJ A«E THE CeEFFKjtATB FOR THE BEarINO CARRIER**! 

C • 

C .. It .ItlJJECJIBN !« ..Mt.MlO# AViP.T)ly:i|(0?4l_ARE.THE CURYE?FIT 

C CBEFp (CIE ntS f6r ThE CBKBI&E6 8t|FFNEBSES IBEaRINO AND CARRIER)* 


IFJ iBSEPBC.CO.CjgB TB 1300 

AKPY. AKP¥«bkPt 7 iaF , 'Y*BKPY) 

AKPZ. ARPZ«6KPr7 (^piieKPt j 

AKTY.AKYY.BKTt/IARTYfBKlY) • .. . 

AKtZ.ARTZ*eKTZ/ J aRtE*BKT7 I 

.. £ 130P_ CONTINUE 

c 

c kPcFx .1 IhpUeS iKvALIo CflWlfSBN " •“ 


r 



OF P<> 




xfchp.? |pfiic8 nFfTP (W/ P8S|*|vi s^Ini 
^PXrB.jjEF.UUS_MF.oi^aW/^tOAUvLlpiNjx«^_ 

ThEPE-U A. DIFFERENCE lN. THE . BAR uCTURE* bF_TW£ BEARING/CAR^I* .. - 
ASSEMBLY SIHUL*Ti«N BETWEEN THf HPBTp AND THE HPFTP 

IF (KpUHP «EG *2 JXp Y'AXPY jKP/ 4 AKPVJKT YIaKTy/xT 2« aKt17kPPY~ , ARPPYJ 

j_KEe4»AKPejiKfclTjAKFJY4KPTl*AttOlZ 

|FJKFO^F.EC.3»KpY»AAPPT|KPZ»AKap:,RTr4AKPTYiKTZ»AKPTZiKppv.AKPv< 

• KFPZ»AKPZ|KPTY M K T YJKPTZ«AKTZ. .* 


«e « AKBC«F.HI0>A«lAKPl*PH|0XA»lPB8> 
— Ct _£_ACBC«FHICZ* • 1 ACB J « PM jDAai ICB 2)_ 
x£A»SKAC4SKAJ.P>,iC)(A*$KAB.PH|0ZS0 
C$A^cA0«8CAl«Pt.|DxA*86Ag.PHjDZSC)-. 

C!A4iCAC48CAl»PH|CxA*$CA2.PHl0Z80 
KSltsKlC4|KU4pH|CxA48Kl2«PH|DZS0 _ 
CSl«4ciC4SCU4PHlCxA455te.PHlOZso) 
_CSij5clCi6Cll«EhjCxAASCZB«EHlDlSQi_ 
K£i.cK2CA8Kai«Pp|CxA4jK*j4PHlOZSG 
CS 2 <SC 8 C 48 CJ 14 PM|OXA 4 S 02 Z 4 PHIDZSO) 

CS£^C2C48CZl#Ph|CxA4SC2B.PHlDZ$o 

THIS SEAfTl N8t IN ?VEt pl^'P 
_j5S3iSK2CA£K3l*EkiCxAJSKJ24pfc!IDl3(L_ 
fCS3«<5Q3C4 8C31.PklCxA4553 e «PHjOZS0r 
CS3 4sC3C4SC3l4Pk|CxA4SC3l.PM|0Z80 _ 


3 CA»(<A0*PH10XA* icAi«PHlDVA4|flA34pn jDxA»bA34pH|DxSQACAA J | )»|t , fAj^V}r<, 
lElEw JX»fl J ».TN8P.i4_PH.lX«EH|X*TNNP| 

IFlFHlX.tTlHThPj) PHIX4FM|*»HTYP| 

SXY 8INIPHXJ ... . • 

CX* C83«PHJX) 

FtPY* PH IDXSQ* ltPPT*Cx»UHPZ*3 X J*PhIDDX« jUHPZ«CX4UHPY*SX ) . 

FLTY. PM|Cx$04(tHTY«Cx«0HTZ*9X*4PHlD0X4iUHTZ«CX4liKTY4S*j 
El8I«_PHl0XSQilUHPY*CX»llHBZ*.SX<4pHlD0X4|UHBZ»cX4UH.et.»S*l 

FCAY. PH|pX3G4(^pAY4CX»UHAZ»SX44PHlDOX«(UHAZ.cX4yMAY»SX) 
FU6lY»PHIDX8C4<liHlY»Cx»WMjZ*SX54p M ip 0x , ( i;hi:, c X4DHiY»S*| . 
Fl.S2Y»PH|0XSG»(ljHZY*Cx.l;M2Z«SX«4PHtDDX« tDH2Z«cX4UM2Y«SX| 

TMS SEAt is h8t In FcEI PUHp 

. ..F«.S3Y«PHlDXSCf.ltH2Y«CX»UMjZ*SX<4PHiD0X«jUh3Z»cX4UH3Y»SXj 

FCPZ* PHICXSC4(GpPY*SX4UKPZ»CXX4PHjPt)X4jUHpz.sx.upPY4CXj 
FUTZ. PHjcxSC»ICPTt«$X4pHTZ»CX«4PH|Ppx4|UHTZ.sX*UHTY»CX j 
Fwez. PhICXSG«(vh b T»SX4U?bZ*CX<4PhIDDX4(DHQZ.$x.uP.bY»C*j 
FCAZ# PHIPXSC«<LhAy«SX4V;HAZ»CxA»PhIDOX4 |UHA2.sX«UHAY»C* ) 

F.CSlMEHIDXSGf IWHlY«SX4UMll?CX;4PHlDOX4tUtnZ#sX*umY*C!j J 

Fl.SZz»PH|P>SC« <pp2Y«SX40P2Z*CX«4PHlDpx»*l/H22. R X.pH2Y»C* ) 


this seal is kst In fuel pump 


A- 


c 

c 



4 


IT’ TV, li 


FCS i | t b»lo • Y • it* VohSI i‘cx5*PH 1 00X« i*UM3 2 ■ *~un3 Y «cx * - 


THIS SECTION DETERMINES 6t*TUS BF SOuEeIe-FILM DAMPING At 

ICAV f C HEAkj NO. SCUEEZC..r?in. DAmpInO 

ICAV * 1 HEAnj SCUEEEE.FJtM PAMNnO WITH CAVITATION ' 

ICASL*_C_hEiNS_SCUEE2E»EJLt_DAHPjNG_kLlHfilH_CA\UJAIlflk. 

f5F?Y»9* 

fsfpj.o. 

F0nY»9* 

FSFTjtO* 

J01CA^t£ftili£liclPJC.8«CiPJL53JN3tiEL6biiELDD2pIEMAiutttAik. 

• iCALl. V 0UB1 lOElbPPtiPELOPPEi 0iypPIi5ECPPl/0*#rSFPY#r8Ffi J. 
•iC1»c1T)C2*C8tJC3*C3TjE1.8D2 , EL^02TjEmAX*ENAXT 

• ice-cut ' 

•ICAIl JOllSl (DECppTY<bELbPTi#Di^PTY/DEtPT2»Of jrSFTY/F8?Tz | 
JFllCA^OtOiCijtl£iC8.CClLLC3j^£jELaDiiEWOD2tlEttA5tlEtlixP_ 

•|CAU.»l8VRM0Et,DPFY»PEL0PP?FDEyPpYtDELP(I:A0» # F3FPY,FBFPi) 

• KiT,cniCclC»T|C*»C3T<EL8D**E>J8DPTJEMAiAEMAxT 

•l£V.*Clt 

• ACALl >)8C/S2(bctcpTY»6EtDPf^Abf u PfYApECPTi»6.'iF6FTY>F§FTfi' 

PgPl4|tEliCELCET*FSFf:. C ILppV ‘ 

^sFT*FSFTY*DeCot^«F§FT2»dECoT9 


OiPfL8*0_£AL£ULAiioN5 

9? I2C9 lPl»lilc_.._ ,-r« - 

EtZe IF(PHlDX,gT»tPHlO>(IPn )08 TO 

08 To *4*1 . 

I F I ?h 1 0 *.» LEI.T PHi 6 ?. Ll P 1 J 100_!8_|2 o6 : 

e6 01 1 1 IF 1*1 

IF(I,EO.b) 00 ToJ2C1 . 

TEFP.lPHlOX«TP^ID>llh/lfpWI6UiPljifPHlb)(tii» 

_ 51.VC»IEP’P« HUY* IP} l*5UYI 1 1 >*5t r ( I » 7 a , i i . 
3L2C.YEMP. ( SLZ ( Jf l )«S|.2I j | }*SlV| 1 1 > ^ 


HaisHisftiai HiwaiiHiKiftiSiti 


!uiS:II": lit*! ! IE ! IzltSJ ! 1 1 l&« ! ! K - *• 


.5%. 


TP 13 SEAt JO n6t HTTuEL PUMP 
|l§?Y0rTErP'IStS3TcIPn*Sl83YlS)) 4 SLs3 Y iI 
SLS320»fEMFMSLS3i(|Pl)«SLS32(8) ) 4 sls32(I 


!?-s*s 




IRIGIMAL PAGE iSI 
F POOR QUALITY 

t 


I 


4 9C. 
-491.- 
4 9c. 

493. 

494. 
<95 . 

494. 
-497.- 

498. 

495. 

7cc. 
7d • 
7C*. 
-2C3.- 
7C4. 
?CS. 
7C4. 
7c7. 
7C8. 
JQS_ 
7}C. 
711. 

7 1 « • 
713. 

’14. 

15- 

7J4. 

?17. 

718. 

7lS. 

7?C. 

• Jtu- 

»?2. 

783. 

784. 

785. 
784. 
282 ... 

788. 

789. 
73C. 
731. 
738. 
-233.- 

734. 

735. 
734. 

737. 

738. 
233.- 
74C. 
741. 
748 • 


TH5 SECTieM CAUCVUAUS $|OElOtDS BASCO.CN. Ak6ULAR-DEvlATt8NS 
8B8H CUPREkT yALVES. T8 fNAK[,E # S ET 161DEL TO NeN-zERfl 


Is 10EL 4 »l ALL SICEL3A09 aDe ZEPO - 

ISlDEt « 1 a L L SICELBaDS aOe RECALCUt* T E6. 


If * ISlDEt •EO*OlQO TO 8150 

SLYT.6LYC 

SLZT.6LZC 

SLYC.lSLYT«Cei<TV*L4THSLejH 

SLZC.«$LZT«C0Htp8l»Th$10S)» 

_SLAt|«SUAtC 

SLAZy *SL A|C 

5lATo»l6lAYT4C0ijTHJL*TVil.9S 

SLAIoil6l.AZTtcesitHSl4TH8l.ei 

SLSlYT'SlSiYO 

SLSltT'SLSlIO 


3LVTtijN|Ti<il4THSLesnt6LK- 
SlPT»3 IN | TH|L*THSl 06 I ) tSLK 


MeLAiTtS|NjTHSL‘*VHSL0sn45LAK 

)«flLAYTtsiN|TM£L4THSL0S)ltSLAK 


-3L8i5oiiSLSmiCBStI^SL#JH3L 

SLSlZ04l9LSllT*C0S{THSL4THSC. 
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B ARRAY 


CC*?? u 
CCOOO L 


R SCAUR 00 *°I l 
B SC»UB CO * 07 L 

. B-. 5 CAU 5 OpC?i ..L. 

B ARRAY OQ 789 L 
B SCAUR CCO* c C 
R SCAUR ccO '8 L 






V 


, . ' . mix ore 

IYpe_cuAss uec, uepps. 


B ARRAY ccttc u 

* ARRAY C 070 A u 

? ARRAY COYeB L 
B ARRAY CC71C L 

.B.>RRaI_ C0Z27-I 

B array 00730 L 


DACj B aRPaT , .00739 u 


B ARRAY 007*2 U 


PPPPPt _ . B SCA 1 .R 0 C 0*>2 u 


B SCAtB OOOCF L 
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-ARRAY.- 007'2 
ARRAY 09O3C 


SCAtR 

$C*LR 

SCAtR 

SCALR 

SCALR 


- 10 - 

10 

13 

15 

-3- 


.!_5CALft_Cc62A-i. 
• " 00*27 ' 

0f«29 

00*27 
00*2* 
0Q7* A 


3 

- 3 
3 

-3 
*0 
— A- 


6CALR OO* 22 L 
SCAtR * 0C7AA i 
ARRAY OO* 80 L 
ARRAY _ 00*^5 L 
ARRAV qq6S5 l 


1 

l- 

1 

1- 

1 


-- ACl 

AC6 

- A!Ri-„ 

BCNOINO 

— CAC* 

CAC0 

— 0AC1 

DAC4 

— OAcf 

-JC4 


ARRAY 00**0 l 
8CALR . 00603 L 
SCAtR 0Q*D5 L 
SCAtR 00*02 L 
ARRAY 00759 L 


> 

.,1 

10 

-10 

10 

-10- 

10 


-1 

20 


I CSi 

-|R8 

!R9< 

VAC l 

JAC* 

-JACJ 


OHMC 

. PHJC 

•PHICBX 

PHjCPT- 

RHJC9U 

_PRJC531L_ 
PHJCT2 
PHJR0- . 
PHjpSf 
PSlCAx .. . 


R ARRAY 006EC 
R ARRAY C070A 
R ARRAY CC76B 
R SCAlR 00000 
_R_aRRaY — 0C72i- 
R ARRAY 0072 a 
R ARRAY . -CC733 
R ARRAY 0073C 
..R ARRAY -007*5 
I SCAlR 00002 
_1_SCAlR_0C*2C- 
1 SCAlR 00*30 
1 SCAlR . 00*26 
I SCAlR 00*28 
I SCAlR - CC7*| 

I SCAlR 007*b 

J_SCAlR OOLAL. 

R apray c07*r 
R 
R 
R 
R 

A 
R 
R 
R 
R 


10 

10 

15 


ARRAY 00078 
ARRAY 0C*3l 
APRaY CC 669 
ARRAY 0C*9F 

.abbay ocfiea 

ARRai 00677 
SCAlR 00*0* 
SCALR 00*0* 
ARRAY cO**i 


L 10 
L - 720 
L 10 

L. ..... 10 

L 10 

.1 10 

L 10 
L - 1 
t 

L 10 


lapel tec 

••••f •• 

ted ccc«e 


LABEL toe 

»•••• • m m m m 

3CCo COCES 


tAett t?c 


LABEL 


. HEX 
toe 


HtX 

label tec 


• ••*4 


LABEL 


HEX 

toe 


local »,ar i able (3 voros i I 

CCCfC 6EKC1 0 CC001 K 


0CC02 I 


LLAKK rOP^OK ( ViCPOS ) 

LAPELEr COUPON BLOCK /BtNCtK/ ji»«3 VORCSj'l 


CCCrC fH|R 

Cc*: 3 -lEE 

CC*?S IRS3 
CCt?F iCSt 
CC**S PhICP 


CC03C PEER 
.CC621.IRI- 


CC**A |CB 
Cc* 3 C ICS] 
Cc **3 PHCPE 


00C78 PHic 

CO 6 2 S !?A 

0C&2B 


ICA 

00(31 PHlCKX 
C066D PHJC*Y 


001*8 PsiCBx 

-0C*2*_JRB 

0062C ICP 
00*38 FhjCBY 
00*77 PniCTE 


00352 PEEC 

0c*27.4?Sl 

Oo*e 0 JCT 
0c**5 PHjCIl 
O 068 I PHJCAy 


00*22 HOOC 

00*28„lRS2 

00*2E ICS | 
00**P PSJCAX . 
00*8B PHICAE 


f 


i® 


CC<*« RHCS V 

Cctcl-fHEE 

Cctr 7 Fhps 
C c7f A ACB. 
CcliA CACS 
C c 7 3 C CACS 
CC7*A -ACS 
CC76S..8. 


CC7*A OFPFs 


0c«9r PHicsti 

cctei-pfciRj 

CC6D8 >Cl 

Cc7i^. AC7 

Cc?2D CACi 
0C?3f OAC5 .. 
Cc7*> jAC* 

Cc74D_AKR 


0C6A5 PHICOJy 

_CC6C3_fH]pJ 

0Q6E2 AC? 
OC7JC CACl .... 
CQ730 CAC7 
OC7*? P*C« 
OC7*C J*C8 
_0C77C_AlJ!l 


boiej Pnici?? o^bo ^hjcbJ 

i)CACi_PUjRB 0O6DS_pHlrtJi. 


00»EC AC3 

-007?1 CAC? 

00733 pACl 

.007*5 0ac7 

007*0 jACt 

_oof»a_A.iB? 


O 06 BD . , 
oo6Ds_pyi 
0oef6 ac* 

0c72* CACl 

Oo?36 0 *C? 

_5c71a_aIr3 


006C7 PMICS3? 

OOID6 .Eh 1RS2 

00700 ACS 
. 00727 CAC A 
00739 0AC3 
... 007*9 JAC2 .. 

007AP OHMC 
00ZA9_flH(!R 


LMFY *0 IMS I 
.CCCfC.fS.APl.fi_. 


f > TERSil. St’EFR ORAKS REOUlREOi 
SjKPLt FEMKC " jPtTVPO 


• ICFES't ERROR EVERITyl 0 tNO 1**0*$) 

CEC. . .. . HEX 

W 0 P D s VBROS 

••• •• • • • • 4 


T8TAL_RR0CRA_t 7*1 flOJEl iELUS_LA?H.IL_CBftPONl. 


> 

OJ 

l£> 


CEAfR/TfC C8C I 76j 

COAfTAM | 1 

LOCAL vAF l ABLE I 3 

TEfP » . c 


CCEr 

coco 

COCO 

ccoc? 


ORIGINAL PA GS ?£ 
OF POOR QUALITY 


CAT* Pf/J.mfiSift!*/ 

CAU. SCUACTDi^/iAC^D.ePB/EiVSCCAl/SALiAl/CoilSOA) 
JFJE.OT.EPAX) E, E PaX 

PE.E.CAL " 

.Fb «E*SAl 

Cl*< l#»E«|)/El.8C2 

Cc »3cPT I t • »PN« • I ) .. ... 

At»A,«(j.«a.je»ci )«cp/i j,*( j»«A.6o«cin 

ASFS • APS (N (PN ) _. ... 

A1*AjaN{(AJ«PN)«|C*FIaA(PK 

flAt*U*;PE/C2J-iAJ»Al.«ASrM 

PS|«At*CAh 

FE'EtCesjoAK) 

PMttSlHGAP) 

CC.1,.PE 

-C1»Cc/EL602__ 

ca»« iso«co«|CRTicC j 

Atlt4?*l?lCl _ .... 

B«.76*PNMI»Aj*4«Cl)rCO 

K*i */Ci/§CRT ( A«i«8«e ) 

JkXiK«C0SlB6U 

Kt.U.SlMPSl ) 

FX»»vSMWX#c6a»VY*56aJ«6T 

Ft»»vS»(WX«8QAfwY»C0A|«Ol 


RETURN 






• Hi f 

* ™ 


*AME.. jtEE.ClA.S LoC WO^QS t«AKI_HPE_ CLASS. 


* ? SCA P occxe 

*«FN R SCA R CCCcE 

L ? $CA P CCCl 3 

cl p sca r cccco 

M R.SCA..F-- CCCcO. 


OEC. ... HEX °EC 

-WORDS toAME IYPE.CLA63 LOC fc5PP«.. 


f 2 

lloc? 

' V 

-olpi 

n 

►SI 

S |K 


CCCC2 L 
OQ0C4 t 
•CCCjO V 
CCOcO p 
CCCd v 
-CCClOJL 

IMP IN 
C<?Oo3 V 
ccoj6 v 

•CCO\A V 


DUMMY 


DUMMY 


.. _•«, 

ATANI 

. CAl ... 

COS 

Tfc 

OAM 

. . Ore .. 

P 1 2 

SCRT 
W _ 

* X 

to 


P SCALR 
R $PRBO 
P SCALP.. 
R SPR80 

_5_8CaLR_ 

R SCALP 
» SCALP. 
R SCALP 
P. SCALP • 
R SCALP 

J*_SCALf? 

P SPR60 
..P .SCALP . 
R SCALP • 
R SCALP • 


QOQ06 V. 1. 

INTPIN 

OOOQ4.V I. 

1&TPIN 

OcOQC-jy L_ 

00003 L 1 

OOOQ* L 1. 

0900 F V 1 

000 SB v DUMMY. 
00005 L 1 

poCCS.v 1 _ 

InTPIn 

00014. V 1. 

0001? V dummy 
OO018. V DUMMY 


ARSiN R SPRoO. (NTRIn „ .. 

A1 P 6CAlR OOOOD V 

COA P SCALP C00U7 V... 

CO R SCALP OOOlj V 

-D1 SLSCAlR OQQOL-I 

E R SCAlR 00002 V 

EX_, P SCAUR «0031C V CUMM 

JOUR* i SCAlR OOOOo V 

PE R SCAlR .00009 v . 

PN R SCAlR COOOA V 

-30A_ ELSCAlR C500B_X 

SOUACT SPRoO EXTERN 

.WX P SCAlP. „00015 V 

XD R SCAlR *00017 v OUMM 


local vapiacl: *23 worCsii 

cccrc *oun. ccoOl pj c6oo§ € 6656} vs OoSoT c^i. 

CCCCE AL OcCO? COA._ ccoca SO* .00003 PE. 0 0 OOA P* 

CCCcC CJ CCCOO At C00CE ASpn O000P Cam Oooio pSI 

CCCl 2 — A CCCiiA ooou « 000 lL_fci 9o01«J<I_ 


00005 sal 
00000 Cl _ 
00011 CO 


tLAKK rOPMeK ( vflRCS | 

L APELEr CepMCN block >6at ijg/' wflPCSiV 


CCCrC Cl 
CCC rE EPA> 


Cocoi C 1 


CCCC2 


00003 C] 


96oo« eLodi 


00005 P I E 


*MRY r 0 IMS t 

CCCfC ~olpT" 


'MR |K<c jC seep OGRAKS USECI 
ARSIN T AK2 -COI 


LxlERNxL EtPFR grams reouireci 


StUCt 


5ATAN2 


•SETUPn 5s in 


5§0rt 


MCPESt ERROR EVERITyI c (NO EPfiePS) 


original pags ;£ 

i OF POOR QUALl.y 



1 

c 

3 

§ 

ft 

i 

7 
£ 
9 
1C 
11 
lc 
J - 


C 

c 

c 


c 


c 


..19 

16 

17 

16 

19 

EC 

cl 
EE 
E 3 
2% 


C 

SL 


__ Iki jo'u* i 1 >b / yd# x # vVo'mb e r n # r y j"” 
JO^flh^J^jfizC Lj Cwp2 i DJiftOD3iiLL£iEtti5L 


.. THIS. I5.TK.?iei9-keEOANCf_rDO.Twt NONi __ 

CAV 1 TAT I K 0 2 Pj PLMN PUll JOCRNAI BEARJNQ 


C£lA_eiZ3i4%18Jf65%/ 

CALL SCUACTU/t.xC^O^OPBiE/VSECALifiALiALiCOAiSQA) 

if « r # ar •Ef.ATc ) .t*znt .... 


EC. !••§?* 

-IB?E2£|?A1j _ 

KE .FpEe.02.cAL/EL/EL/S0Rt (EC< 

V0fppO£.8AL/EL/SC«T(ECl 

.>V*E!CAI"*B.5aC-.- 

WY.kE«SAl«WB.Cal 

f*..v3« (**«coa.*y.$oA )#bi ! 

Ft*.V9MHA«S0A4^T#C0A)#01 

RETURN 

END 



930 


932 

933 

93% 

935 , 

936 

93Z 


9%0 

--*11-.. 

9%2 


9%7 

9%« 

9%9 

950 

961 



y> 

-u 


P>PE — i»PE -CLA-S- 


HEX 

-LeC- 




>L 

Cl 

C 3 

lc 


r sc* 
P SCA 
P SCa 
P SCA 


CCCC4 V 
occcc L 
CCCc? L 

CCCo A V 


Lpax R. .SCA- R — CCCc&.L. 


I Y 
tPP 
Jai 

SCLACI 

h E 


SCa 

SCA 

SCA 


p *ccci7 v duphy 
R «CCC|S V OOP.r.Y. 
R OOOcS V 1 

. SFR a ..EXTERN-: 

P SCa p CCCcC V t 

B -SC A_ fl_«0C0i3_y_0L'Rh*. 
P SCa R *CCCi2 V CUPPY 


-RAPE — JYRt— C laSB- 


HEX 

rtflC- 


OEC. 

-JOHCg. 


CAL .. . - R SCALP- C0OOA-V i COA _ 

PI R SCALP OCCOI L 1 0? 

| * R SCALR--OOC02 V 1_ E9- fs 

IL R (CALR OqOCO v i EUOb* 

-Ea R— $CAL b — C0CC9-Y 1 PX 

J8U«* 9PR68 00030 p jO'jftt 

Pi R SCAUR -0OOO1 -V l-_ PI*—. 

SO A R SCAUR ooooe V I SQR'f 

VS ...... R. SCAUR — 00003 V I — WB 

*X R SCAUR OOCOf V 1 WY 

JO R.SCALft— »OcOll -V -t)UP.fc* » 


AJAmC — TY pE-CtAsS 

. — . . 4 • . — — . . . . 

R SCAtR 
R SCA(_R 
R SCAtR 
R SCALR 


HlX 

-U8C- 


OEC 

-kBRDS. 


00007 V 1 

00002 L I 

OOOOQ V 1 

ooooa u I 

-R— SCA^R aQQOl 4— V— CUPMY. 

I SCAtR OCOOO V 1 

R SCAUR 00005 U 1 

R SPROO JNTHIn 

-R SCAUR -OOOOE V 1 

R SCAUR OOOJo V I 

_R_ SCAUR-* 000 lk— V- 0 UPJ 1 T- 


LOCAU wiPlAfiLE (17 keRDSlI . . 

CCCcC_wBUP2 DOOCLJ!] OOOCA-I 66§03_V$ C 660 S_CAL Y)000a_SAL 

CCCr S AU 00007 COA C00C8 SOA OQOC 9 t? COOOA £C 00000 Eo 

cccrc *e ccooo iu „ . Ooooe *6 ooccr >x Oooio wy._- 


PL.'M rBPEflK ( W6PCS | 


LaPEuE? C8PP8N BuOCX /OATAjB/ (7 kOflCSM 

cccrc ci. ooooi oi 666c* 6* oooos 6j . Ooofr clb6I" “ " 666ob p I* 

CCCft ErAX : ... 


emry abimsi 

CCCrC *8U'P2 

iMRikiiC SLeP .eGRAPS .U3ECI 

scr’i ...... 

exurmi JteFR or aps recuIre'ci 
SCL»cT ~ _ $etlfk TscrT - 

h IOPESt ERRS 9 EvCRJTyI C (A8 EPb; p S j 

CEC HEX 

kCRCs V 0 RDS 

*•*?* •••?? 

CEPERaTEC C60 i 95 COO«F 










\ 


MME-..»YtL.CUA S. 


Mn , -TTEt-CiUS-^fic^-tteROs MHt_^iy#c_cUAss Sc hbpds 


C3 

IMA* 

*Y 

► 

-burs 

*12 ..„ 
S|K 


B SCA R 
R SFR 0 
P SCA R 
R SFR 0 
.£ SCA_B. 
R SCa A 
* SCA A 
R SCa A 
R SCa R 
I SCa r 
-R.ISa.B. 
R ifr o 
R SCA R 
R SCa A 


CCCC? V l. 

IKTB |N 

CCCo*. V. ... I 

IMR|N 

_cccct_y l 

CCCca U 1 

CCOCS l l 

•CCCjO V DUWU 

eccis » t 

CCOcO v 1 


. R SCAUP-OOOU.V. 

R SC*tR ocooe v 
....P SCAUP _OOOOA V. 
R SCAUP 00002 V 

fL.SC ALB OQOCl-l- 

R SCAUR 0C005 V 
-„P SCAUB_ O00i0_v. 
p SCAUR 00009 V 


Alt 

BO 

CU 

Cl 


P SCAUP-00012 V 


_.Cfi9o5_V A. 

!Mr|N 

•ccci9. v oumvr 

•CCCjA V DUMMY 


R SCAUR *00016 V burn, 

JLiCAUR ooQCf v . [ 

R SPR60 INT R IN 
. P . SCAUP . *00017_V DUMMY. 
R SCAUR *00016 v DUMMY 


It 


fl8D* 

a— 

mn_ 

HO 

xo 

YO 


_ P SCAJ.R 00018 V .. 1 

R SCAlR 0000% V 1 

._B. SCAUR 00000 U-. I 

R SCAlR 00003 V 1 

B_SCAlR OOQU? t 1, 

R SCALR 0000% L 1 

_.P_SCAlR *0001C V dummy 

R SCALR 0001} v 1 

SPROO . OCOOO P 

R SCAlR 00001 V l 

R S CALP OOOOB V 1. 

R SCALR OOOU V 1 

_B.. SCALP ..OOOOC..V 1 

R SCAlR OCOOO V 1 


c e C AL_viRjA?lt t SUiOROSj 


ccccc WBURJ 

CCCf 6 c 
CCCtC >C 
CCCtS ► £ 


CCC01 Pt ... 
0CC07 A 
COOOO YO ._ 
CcOll M 


Oooos Cl 

00CC8 h 

OOCQC 0 ... 
0001% VO 


ooco? c? .... 

00909 OAMO 
OOOOP 0. .. 
00016 AK 


?Sooi cSaR 


$0010. to 

Pool* At 


00005 c 

OOOOB SOAH 

....00011 00 . .. 


' LAM fOMPOR I K8RCS ) 


CAftLtn COMBM buocr /DatajbA . it. hBrCSi I 


CCCtC-CU 

CCCf* IMA* 


_C0091_D1 


i000t_O6. 


5&0l_0i. 


-OoQ&li^BPt DOOOajLll 


•A!Ry JelMSI 


CcCf.C..vOCRj_. 


iMRiMic suep borams uitoi 


• COAT 


MttMiL SueFR ORAMS BtOU|RE0l 


SaTaM 


9PVRRR 


95ETUPN 


•SIN _ 9SCRT 


►icMS» tMofi'"tvlRitti C (KB tRfi8*sr 




<3 


tit ’ 


ORIGINAL PAGE 19 
OF POOR QUALITY 



m 







i 


j*fE_£lA-S- 


HE* 

-LOC. 


n F SCA B OCC03 

Cal ft SCA ft «CCCoD 

cl f sca r cccco 

c; ft sca ft cccca 

CpAX S-£CA_ll_CCOcA- 

S *L R SCA R *CCCcE V DUMMY 

SCftl R SFft c I k t r i n 

•*s R SCA ft *CCCcC v CLMHY 

'CS ft SCA ft .cccol v_. 

>C« ft SCA R COCca v 


CSC. 
-keftog. 
• •••« 

.... I 

DLMMY 

-- I 

I 

u 


HCx 


oec 


1 



r scAtft-«6pocr. 

ft 6CAlft 0Q005 
R -SC*Lft— OOOOl. 
ft 5CALR 40C00B 
ft-SCALft— aQqOOA. 
ft SCAlft 0600 % 
ft SCALft - OCCOO 
ft SCAlft «0$000 


V- DUMMY. 

V l 

L 1- 

V DUMMY 
V— OUHtlt- 


, . OEC 

*AMC— Urpt-CUAgS LOC N0RD1. 


V i 

V 1- 

V DUMMY 

.Y _ ft_SCAUft-«00007_v-0UhHY- 


-ATAN* R SPftOO -|NTR|N 

COA R SCAiR .00010 V CUMMY 

Ct 5 R SCALR 00002 L 1 

EL068 R SCALR COOOA l 1 

-ft It R—SCAift — 00008— I 1- 


90 A R SCALR *00011 v DUMMY 

60UACT SPROO -OOOOO ft- 

XO ft SCAlft .00008 V dummy 

.YO R -SCAiR .00009 V OUMMY 


NEx 

CAEti .coc 

••••f % ••••• 

i' cccac 


.. HEX - 1 , HEx : — -Wt* - 

label ioc tAeti ioc label loc. label . loc 

••••% ...JM gtgftg • •••g % f ^ % 

a coon. 


LABEL lbc 


HEX 


• £ • ••••£ 


cn 

o 


local variable <• vbbosm. ... 

CCCfC SCLAC OOOOl X09 


ooooa .T?9.. 66 oo3.ae .• Ooo6i itoET oooos cbet 


FLAKk rOPFBK ( NORCO ) 

LABElEr COMMON block /DMAje/jl' WORDS fl 

cccol ci . oocoa ca 


cccfc cl 

_CCCf.t_Et>X_ 


00003 C3 


(LOOS 


00005 P|t 


entry foimsi 
CCC fC SCLAC 


INTRINSIC SLep OQRAMS USEDl 
ATAnc CRT 


2 ° 

53 

TJ O 

o ~ 
o 

» r 

O -s 

-C > 
> I s 




External sieRfi gramr reolireoi 

SaTaNB SETLftN 9SCRT 


► 16hESt-£PR0P_EYEP.!TYJ..C-«N0_ERftOftSl 


CEC. 

WORDS 


...HEX. 

WORDS 




assign (Mooio /recw'i 

t*Le*a_;oaJajE4-£4££aLtJu4uoc|A4iJ^ni»t»iooaiJU4i*iti*i» 

IllB/i SE«#$YS Er)M^AP»ALU 

*SS ION iFllie, 2/DM1) 

*SS ION (M 1 It/ 2, CAT?) 

ASSIGN IFIICI/ YACl)/vrC 

assign (Mice# yaci>*vfc 

ASSUN-4r-U22/-J4CU/-vEC 

ASSIGN ir«5/Sl to) 

ASSIGN IFI2CC/ 2 #eENO» 

ASSIGN |F I |2 1/ Fa02) 

ass Ion . if I3CI/ t/poct) 

ASSIGN (FI3CE/ 2/TOC2) 

«SS IGN— (EJ3C2/_2/C0C3| 

ASSIGN (F>3C«/ 2/PODA) 

ASSIGN IFI3CS/ 2,)*0C9) 

ASSIGN |F«3C«/ E,rOC6) 


ASSIGN CFI3C7/ 2/P007 i 






h oeR#f p-A* 

''PEORO.'kC.-FJOQ. 


nr fF#fiE)i 

K uee* or eec* ms« 1 root., . 

c cophOn, 

{L.OVERUx 

. . c public LjepARY 

MMTSi FWAI ECCC |>AI 7007 _ 

L'etKC « 

H l*M> CO^OR.. ASE SUE . 0 

r Letic LieaARi s kone 

*CIAL..^EIt 6RI.8_2EJ 5008/ 0301 0 , »fl RDi 

'DIAL PILE £12 I r*/ (DC QRakU.EO 


UlAL.^EtBRI.OJEJ £008/ 0301 0 , fcfl RPi 

' 0 1 AL PILE £12 I r*/ Oec ORpkUlEO 

nn& ias< -imin-T 


IVIOAC Fte ENTRY UISPhA SUE BtSFHA SIVE DCOfiB 
7<(i . Eccc OCCE _ 779 a 3E»..._ 7B»0 96__.77|A 


■ 0T_7il 


iKFLir UMM? Wp-fii. lP"lt" 7 SX“'~lW '"~tWf'^ 5 HiRS"LwSR‘ ' 6 ffXN' 

APS? 13?C S0C8 . lfp<8 IOCO.„ ?BOf t*CE 6. 0 t 


CCMFkL SSCTIM. 


























ORIGINS- * 

OF POOR QOAL' 






















/ 


u 


i * 


cn 

00 


Ik 


err u 
-ctr.-tt 


ccr u 
ccr it 
ccr u 
ccr u 
ccr ll 
ccr.LL 



• Inpott 
cr.soer-f. 


cioocerc 
7ST0P .... 
9 1 TOO 
90TOt .... 
JStTUpp 
7StI_ 


COO 

ippt*cr 

|RTSSn« 

ITO 

SRCxh 
-7PHC- 


7CRRIMT 

tOCBNApC 

|abb"teA 

IWRCCxjT. 

tABRTiev 

KBCExjL 

ICRRExtT 

criBVTRO 

CtOTRiO 

CSCNtnC 

CTINlT 

|A6vArtB- 

70CBSRCH 

(OBDCiaC 


If PTC mn 
irpjPiR 
JREAtL 


7terA8RT 
IVINITvaI 
C TLB 1 T 
JOCOlt,.. 


IKRCtPTl 
7 JBOAIVI- 


leoiive 
IIOOAcoR 
CIOOTtpC 
IIBCTa 
SNBuTcjV 

tcholoi 

lINtJcpt. 
neiNtoc 
itcoIt . 

iTRKCTf-F 

ITRUNr 

scata 

_*OAJI|kR.. 

7 ERRBr 

..IOMTkaP 

70NITAC 


T _ 


o o 

"n -X 

o ^ 

o > 

S3 r* 


*CJ 


A 


•— m **in tw 


CIF U 132 C 

CEF-.U. C3E.C. 

CEF LL 3C7 C 

CEF LL 22? C 

CEF LL 2)5 C 

CEF LL FFFF FPQ C 

CEF LL 2CF C 

CEE .LL Z_C_ 

CEF LL 1EA C 

CEF LL 1E2 C 

CE* LL 2 1 A C 

CEF <.L 25? C 

CEF LL 259 C 


CEL LL 23J_£_ 

CEF LL 25A C 

CEF LL 3C3 0 . 

CCF LL 3C* C 

CEF LL JC? C. 

CIF LL 3C* C 

Cel. ll 

CEF ll 3C® C 

CEF LL JC? C 

CEF LL 3CA C 

CEF LL 3C9 5 

CEF LL FEE C 

CEE-LL 31*_C 

CEF LL cC£ C 

CEF LL 2E9 C 

CEF LL 3jC C 

CEF LL 
CEF LL 


CEL.LL 351-C 


CEF LL 
CEF LL 
CEF LL 
CEF LL 
CEF LL 
CEF_LL 
CEF LL 
CEF LL 
CEF LL 
CEF LL 
CF.F LL 


ME C 

51* C 
MA C 
??E 9 
8 fJ 0 
5E6_fl_ 
88* C 
(89 C 
63E C 
9 ?f a 
8(9 C 


C*ClNC..k*I C8t_LEJEC 

TILE FMFLEX lseo ibo.oaakulcI 

FILE 2t*Xl ueeo . e g*akvleb 

file ._eT#x? t$EO 4_6"AMJLE$_ 


JUNOEC 

.IHtblER 

•EOFt . 
HRBflRt- -- 
8REKTBLX 
IPIB9FF.. 
IVAnSipT 

-IbBABwci 

7BOF8LT 

JBOFOutC... 

8TEA«8R 

• iFCykO. .. 
•TDI*Ig1 

-UOJElfl 

86|nBof 
.(EBFCalI 

CEOFJacR 

• EOFOac" .. 
ccorutn® 

JPECfU? 

(WreCalI 
EWPEJACg ... 

«kreOac r 

IOFTIrp 

JlCVjE_ 

•tiNPcrt 

UPNCjpl 

•CVR^T 

j^bet. 

CIS . 

J:EPJ8|I 

RCkD 

CIBt r 

KI9F 

.KPA,...- . 
ERRACp 

-abuaBo 

ERTAO 

CBHSTr 

;«h 

DISCBk. 

enocBk 


FILE 5T»>3 



3 O^ANVLES 



or;: 

OF POOR Q'JA 











ics Jicr 3/7 .Rer #PERAT|Na.tY«Un .. ?• ..YlRStON.Coi- ••. _ 


tatH »e t I f'C • ••• ociitiet 





ORIGINAL 
OF POOR 






>cs sjom •✓? bbh BPtfuiifcj System VTHsteN. c6i,_ ••_ 


n"i •(••'•'••• 

aCC®U^T •!••••!••• 

ieT>L si TIME OCfOiTJ? 








umuinMi. raot la 

OF POOR QUALITY 


: 12 : 6:03 

TIME 10 

• T <1 

* w ) 


n y vnui a a 

0000000 

0.0000000 


0 . 0000000 

r. . 

0000000 

0 . 0000000 


0.0000000 

0 . 0000'' 0'' 

1 100000.1*02 

* 9 • ^ * # # * .« 1 « 

~r 4 

1 4£ 400002-05 

1 , !!• 1 Z ? 9 ? 

0000000 

0.0000000 


0 . 0000000 

0 . 0 0 0 0 1* 0 o 

o nv??w^£~ o ^ 

3.1 500000E 

— a r 

w w» 

2 . 2729 9991" 0" 

1 .3910000 

ooooceo” 

0.00 .'0000 


0. 0000000 ' 

0.0000000 

0000000 

0.0000000 


C. 0000000 

0 . 0000000 

0000000 

0.0000000 


0.0000000 

0 . 0000000 

2 0 0 0 0 0 0 E - 0 3 

1 ♦ 5760000E- 

-03 

1 • 09700002-03 

5 . 4 330010! 


2ZT-: - 4 .9999999E-03 

~ -30505 » 23 O -Wi nonrwn&l 

roots: 

-0 . 7993591 1E + C2 0 . 44 ?2 » 3 ‘,*£f 0 •’ 

■ 0,77? 35 HE + 02 -0 .44724367E+04 
■0.27631234E+03 0 . 4 0 9 6 1 4 9 2 E + 0 4 

-0*2762 1231E+0 2 -0.40961492E+04 
“0.860663S1E+02 0 . 4171 7931E+04 

-0.H6066981E+02 -0 . 4 1 71 79 3 1 E +04 


"0 * 7352‘234>3E + 02 
“0*7Ji28846E+02 
■0*4c71 6722E+03 
-0 . 4351 6722E + 0 i 
-0.23956336S+02 


0 . 4104341 4E+04 
“0*41043414E+04 
0.34113553E+04 
-0 . 34 1 1 3332E+04 
0.37 4 9 2 9 3 9 E 0 4 


33*563363+02 

-0.234929393+04 

289319301+02 

0 . 30 *30 377E+04 

25 95 1 9S0E+02 

- 0 * 3 0 4 5 S 8 7 7 E + 0 1 

15S83320E+02 

0 . 253963 1 1E + 04 

19S83520E+02 

-0.252963113 + 0 ! 

21216342E+02 

0.272648153+04 


0 . 21 2 1 6342E+02 -0 . 2726481SE + 04 
-0.23704754E+02 0 . 52636377E+03 

0.257047543+02 -0.52636377E+03 
0.40716235E+01 0. 687762 75 E+ 03 

0.40716235E+01 -0 . 637 76 275E +02 
0 . 1228°674E + 02 0 . 2736623CE+03 

0.13289674E+02 -0 . 27566230E+03 
0.10169385E+02 0. 22225 S77E+ 04 

C. 10 16933 5E +02 -0. 2222337 71 +04 
0 . E3239346E+02 0. 16278469E+04 

0 *83239i^6E + 02 -0 . 162 734 6 9i£ + 04 
0 . 56057606E+01 0 , 20655027E + 04 

0 .56087606E + 01 -0 . 206330 37E + 04 
0 . 2 1 1 2 4538E + 02 0.1527 6960E+04 

0.21124538E+02 -0 . 1 3276960E+0 4 
0 . 5 7=41 4 6 1E+01 0.19201467E+C4 

0 - 9 4 1 46 1E + 0 1 -0.193 0 1 4 £ 7 Z + 0 4 




«TE 21 9184 TIME 16J13 

! t' ♦ 

• Ilf 

0.0000 000 0 >0900000 

0.0000000 c .0000000 

2.1100C00E-02 9.434°??l£-04 

; ^ . 

» «# # 

0.0000000 0.0000000 

9.0279797E-04 3 . 15000001-02 


ORfGfNAL PACt! {•* 
0F POOR QUALITY 


c . 0000000 


0.0000000 


0 . C 0 0 0 ' 0 0 

0.000000c 

1.46 400002-05 

0 . 0000000 

2 . 2‘>S7?99E-03 
0.0000000 


a r r r ijiinr 

0.0000000 


0 . 0000000 
1 .8710090E-03 
0.0000909 


0.0000000 o.oocoooo 

0.0000000 0.0000000 

2 . 2 0 0 0 0 0 0 E - 0 3 1 . 57600002-03 


0 . 0000000 
0.0000000 
1 .0970000E-03 


0.0000000 
0.0000000 
9 . 4530910E-04 


ZER = 4.9999999R-03 

AK : 

2750000. 2750000. 

-5*5.0540 -535.0540 

1 . E237799E-02 1 . S237799E-02 -<! 

9.2844002E-06 9 . 284 4092E-06 : 

OM « -30505.23 6A0= 0.261! 

racist 

-0. 1 1840662E+03 0 . 46 422549E + 04 

-9.1134C662E+03 -0.464223492+04 
-0 . 640307S3E + 02 0 . 4 1 S 3 ~ 6 5 5 E f 9 4 

-0.64050783E+02 -0 . 4 135' , £53E+C4 
-9*54038 3 3CE+92 0.41118126E+04 

-0.649S8330E+02 -0*4’ 11812&E+0 \ 
-0.23978426S+03 0 . 39 1 73275E+04 

-0.25973426E+03 -0 * 371722752 +0 4 
-0 . 43527661 E + 03 0 . 34 1 1 1217E f O 1 

-0. 43527661E + 03 -0 . 34 1 1 1 2 i 7E + 0 » 
-0. 38991 S96E+02 0 . 333477932+04 

-0.3877139 6E+02 -0.333477932+04 
-0. 29878 0S5E+02 0 . 304065961+04 

-0.29973085E+02 -0.30406576E+04 
-0.23020304E+02 0.29192225“ +04 

-0.23020304E+02 -0 . 2? 1 922832+0* 
-0.22762785E+02 0 . 27265006 2 +04 

-0.22762735E + 02 -0 . 27263006E+0* 
-0.38109047E+01 0.63C24130E +02 

-0.33107047E+01 -0 . 68S2 4 1 30G+03 
-0.256 4288 1E+02 0 . 52692359E+0 3 

-0.256 42881E+02 -0*52692389E+0 3 
-0. 132S3338E+02 0 . 27366540Z+03 

-0. 132S333SE+02 -0.275665402+03 
-0. 1 1775142E + 02 0 . 22 1 75 l ’9?E + 0 4 

-0. 11775142E+02 -0 . 22 1 " , 597?!I +0 4 


4050000. 
-814.8732 
6 .22222E3E-0! 
2 . 9 0 2 6 4 0 0 E - 0 ! 


4050000. 
-814.8732 
■ i .33322332- 
2.9026 4 COE- 


; 733 i - o : 


; T - 1 T O T *. ( 


-0 . 4730230 1E+ 02 
-9.4°392391£+92 

n 710 flOOj''Clrt^ 

/ I V. “ . w T _ •_ 1 V w 


-Or !. 1 

C , 184752 * 02 - 0 * 

-A i 2 4 ** f J '• ') A » 

M * i vv ..•> *+ — . • 


C • 3 430G742E+00 -0 . 17608734E 
•0.44176924E + 01 0 , 192221232 + 9* 

•0.44176924S+01 -0 ♦ 17222123E +0 • 


B3 


r« to 
O O 





ORIGINAL PAGE U5 
OF POOR QUALITY 


# ,T? A • O • O A 

4. • 

TIME 15: 55 


*• • • # 

c . oococoo 

A fk A A A A A /) A f A ^ ^ A A A A 

A # A A A A A A *5 

C. 0000000 

<\ A * A A A 0 A A A 

c . 0000000 

2. 1 1C000CE-02 

9.43499913-04 1 . 46 400COS-05 

i ^ i joo.vir 

er* • 

h 4 * 



o.ooooeoc 

0.0000000 0.0000000 

0.0000000 

9.0299939E-04 

3. 1500000E-03 2 . 2799999E-03 

1 .E910000E 

0,0000000 

O.OCCOCOO 0.0000000 

0.0000000 

c n • 

mU 1 



0.0000000 

0.0000000 0. oococoo 

0,0000000 

0.0000000 

0.0000000 0 .0000000 

0.0000000 

4 .5100000E-04 

4 .9999999 *3 - 03 1 . 0970 OCOE-C 3 

9 . 453001CE 

E-ETA = 0.0000000 


ZER * 4.9999999 

E-03 


fi k: 



2750000. 

2750000. 4050000. 

4050000. 

-525.0540 

-535.0540 -81. .3732 

-314.0732 

1 »S2377??E-C2 

1 ♦ 223?7??E-C2 -6.3SZ222ZE-02 

-6 . 32322232 

9 , 23 14002E -06 

9 . 23440023-06 2.90264003-05 

2.9026400E 

CM * -30505.23 

Ga£ nonV*a&l 


FOOTS : 



-0.71952374E+01 

0 .447399453+04 


-0.719672743+01 

- A ^ + “7 o o o {j »s; r jl r. 


-0 . U978194E+03 

0*41001 4 70S + 0 4 


-0 . 439781945+03 

-0 >41 001**70 EvO*?. 


0 . £?2S04i4E+00 

p , 4 1 iosA*:crAA 4 i 


C. 692204142+00 

— a ^ ^ 13203 9 3 10 ' 


-0.63169950E+02 

0 .410395063 *04 


-0.63169950E+02 

-0 .410575063 *04 


-0 ♦ 485281 30E+03 

0.34113952E+04 


-0 . 425221293+03 

-C.34113953E+04 


-0.33347678E+02 

0.325137543+04 


-C.73347673E+02 

-Q.33513754E+04 

• 

-0. 297524 03E+02 

0. 304602 01E+0 4 


-0.297524033+02 

-0.3046020 IE +04 


-0 . 21 ?830?2E*f 02 

0.29Z95590E+04 


-0.21933092E+02 

-0.29293^703+04 


-0 . 225 1 1638E+02 

0.27259303E+C4 


**0 » 2 23 1 1 6 38E +0 2 

-0.2-2693033*0. 


-0.13281420E+02 

0.275661373 103 


-0. 13231 420E+02 

-0. 275661373+03 


-0 • 2363S004E+02 

0.52634 6 1CE+03 


-0. 256250045+02 

-0. 5" > * °4 Aiupifti 


-0.3S795707E+01 

0.62753533E+C2 


-0.33795707E+01 

-0.62768w53E^03 


-0. 11221522E+02 

0 . 222Z0394E+04 


_r, _ i io*»i =• t ">r- 4-0 

-0 . 222202° 4ZfO- rt 


-0 . 660621 403+02 

0 ♦ 1 1 2 5 ° i -i 0 Z 0 4 


-0. 660621 40E+02 

• 0 » 1 .?• 25 ° ' l ? E 0 


-9.11334725E+0Z 

0 . 20 >7 c ?4^2E + e 4 


- 0 . 1 1 2 3 4 7 2 A 7 - 0 2 

— a ^ a f . j ^ r rr c a 



0 > • 220? •*:;■£ W'4 


-0 > 20224c77Ef 02 

-0 . 122^3 ’^175 


- .> 1 A T “7 O L t "\ T - /*i ? 

a , 1 ° 2 C 1 0 0 2 E ^ 0 4 


... A lATTO.iTfMTiO 1 

. * It / W -J W — 

- A lO*JAt AA'CLA 



t J ro o O O O O O O 

M >1 o o o o o o o 

10 W O MOO o o o 

w rj ♦ o o o o o o 


re 


2 : ?:s4 


TIME 


original n " r " T<7 

OF POOR Q.'-A' '17 


k: 

0.0000000 

0.0000000 

0.0000000 

0.0000000 

0.0000000 

0.0000000 

0*0000000 

0 . 00000 ‘ 0 

2. 1 100000E-02 

, * 1 

9.43490915-04 

1 . 46400005-05 

1 . 2159999E-02 

* 

0.0000000 

0.0000000 

0.0000000 

0.0900000 

9.0297989E-04 

3. 15090005-03 

2.2739999E-02 

1 .S910000E-03 

0.0000000 

n: 

0.0000000 

0.0000000 

0.0000000 

C. 0000000 

0.0000000 

0 .0000000 

0.0000000 

0.0000000 

0.0000000 

0 .ooooo'oc 

0*0000000 

4 . 5 1 0 0 0 0 0 E - 0 4 

4 , 9999999E-02 

1 .C97C000E-02 

9.4520010E-04 

’ETA = 0.000000 

0 



ER = 4.9999999 

KJ 

5-03 

« 



r> 7 n rs r> r, 

-535.0540 
1 .8237777E-02 


9.2S44002E-06 
0M = -30505.22 

F.QCTS : 

0 . 45324308E+01 


4050000 . 
-814.3732 


2750000. 

-535.0540 
1 .8227797E-02 -6.3S22233E- 
9 . 2344092E-06 2.9026400E- 

C-A<3 = .l<el5050 


02 


05 


4050000. 

-314.8732 

■6.3S32283E-C2 


2.9026400E-05 


0j45 4 9 31? 0 E + 04 


0 . 453243035+01 

- 0 . 464731305+0 4 

- 0 . 321381635+02 

0 . 417891805+04 

- 0 . 3213316 GE +02 

-0 # 41 7 G? 1 8 CE + 04 

— 0 * f) 212 £.*^ ocr 'f 0 r 

A . .44 «• 

- 0 , 521266435+02 

- 0*41 1522 :J 2 E *-0 4 

- 0 . 419063515+03 

0 . 371801445+04 

- 0 . 41706351 E +07 

- 0 . 37130144 E +04 

- 0 . 435905295+03 

0 ♦ 3 4- 1 0 9 0 C 5 E + 0 4 

- 0 . 4 S 390529 E +03 

- 0 » 34 1 C 7003 E +04 

- 0 * 425217073+02 

0 . 33400184 E +04 

- 0 * 425217073+02 

- 0 . 334001345+04 

- 0 * 305325533+02 

0 . 30410829 E + 04 

- 0 * 305825533+02 

- 0 . 30410327 E +04 

- 0 . 24395551 E + 0 2 

0 .292003685 + 04 

- 0 * 245955513+02 

-C . 272003685+04 

- 0 . 2394902 SE +02 

0 . 27274074 E + C 4 

- 0 , 239490285+02 

- 0 . 2727407 4 E +0 4 

- 0 ♦ 23587810 E +02 

0 . 52690550 E +03 

- 0 * 255078103+02 

- 0 . 52670550 E +03 

- 0 * 349755473+01 

0 . 68817632 E +03 

- 0 . 349755473+01 

- 0 . 63317632 ET +03 

- 0 . 132759793+02 

0 . 273662295+03 

- 0 . 132759793+02 

- 0 . 273662275+03 

- 0 * 120953133+02 

0 . 22177543 E +04 

_ A innoc-riTCAAT 

w > A 4 « w / W W 4 w L • w .« 

- 0 . 221775 ^ 32+^4 

- 0 * 437693 SSS +02 

0 . 172756255+04 

- 0 . 43769333 E +02 

- A 1 1 ' 1 -iu < n b- c A A -1 

y • * / ■ m * W . 't 

• 0 . 41321 1 3 2 r 4 - 0 ^ 

0 . 13538065 E + 0 4 

- a r a 1 2 2 1 1 3 2 £ 0 2 

— a , < o g . % ’.‘■•rTA.y.j 

“ 0 . 1 102491 5£ ^02 

0 . 1 96337 • 9 £+ '■& 

-C • 1 i 024 o 13 £t 02 

- 0 . 176337155+04 

— a # 97 i 322 n 7 EiC'l 

0 . l a 29036 '>£ + 04 

-a ^ ?^£ 2322 7 E i *0 1 

- A mi/NAc • r* i. /\ 4 

V. V to 0 — * m mt -J w *... 1 .* . 













ORIGINAL TAGS iS 
OF POOR QUALITY 


. 


I. 


» • 
1 . 


• • 


TE 2J 91S4 

t.- • 

TIKE 1 £ 125 



0.0000000 

c 0 a a a r, r A 

A # A A A A A Q A 

0 . A A A A A A A 

0.0000000 

0.0000000 

0. 0000000 

0.0000000 

2.1100000002 
A « 

9 • 434???1I~04 

1 . 46400003-03 

1.2139999 

w • 

0*0000000 

A r AAA AAA A 

0.0000000 

A A A A A A A A 

V f W »• • to » w 

9.0297929E-04 

3. 150000 CE-03 

2.27399993-03 

1 .3*10000 

0 . 0000000 

n • 

0 . 0000000 

0.0000000 

0*0000200 

wt • 

0.0000000 

0.0000000 

0.0000000 

0.0000000 

0*0000000 

0.0000000 

0.0000000 

0.0000099 

4.51000006-04 

1 .5760000E-03 

7.10000003-03 

9.4530010 

ETA « 0.000000 

0 



3?: = 4.9999999 

W* • 

E-03 



. » 1 

273C000. 

2750000 . 

4030000. 

4030000. 

-585*0540 

-585.0540 

-214 .2732 

• -314 .3732 

1 . S237799E-02 

1 .62377993-02 

-6. 38322836-02 

-6 . 3E322E3 

9 . 28440026-06 

9 .28440026-06 

2.90264003-03 

2.9026400 

M = -30505.23 

GA0~ Q . 2 i 

15050 



0. 113-93564E + 03 
0.1 159236 43 + 03. 
C . 99 f 43399E-C 1 
0.9?5-»32lT?I-j'f 
0 . 1219 C 730E-K3 
Ctl21??7 , 0Et33 
0.262124963+03 
0. 2621249 6E +03 
C.43527930E+03 
0*485279506+93 
0.39719165E+02 
0*397191656+02 
0.297737663+02 
0.297737663+02 
0.23175104E+02 
0.231751043+02 
-0.25391 1633+02 
0.23391165E+02 
0.371154193+01 
-0 .371154193+01 
-0.25593602E+02 
-0.255936023+02 
-0. 132760526+02 
-0 . 13276032E i-02 
-0. 119812846+02 
-0.119S1234E+02 

A > 4O5O4A20E + 0'* 
0 . 4 ? ? 9 £ 0 0 ” + 0 2 
0 . 47747221 E + 02 

a > 4 ' 7-?4 — o 15 1 ^ 0 ° 

0 . 4 — 7 — •230Ef'*' 1 
a , 4272 4 230E +0 1 


0 . 463S4316E+04 
-0 .46334216E+04 


C • 4 1 £ 0 2 0 3 7 E + 0 - 

:• 41 £ 02 V.-j 72 •••>-: 

0 . -ri361221E + C J 
-a . 417612212 !- 0 4 
0 . 292232373 -C - 1 
-C . 39225237E +04 


0.341112133+04 
- 0 .341112152 +04 
C . 333312C3Er04 
-0*333512036+04 
0*304066456+04 
-0*304066456+04 
0.291°3743£+04 
-0.29193743E+0 4 
0. 272741243+04 
•0*27274 1246904 
0.6S817047E+03 
-0. 63317 0 47E + C3 
0.52690164E+C3 
-0 .526901646+03 
0.27366200E+03 
-0 .273662003 +03 
0*221708556404 
-0.221768556404 
0 - 17333^963404 

— A 1 7 1 t v \ r a , 
r # 1 *3 ^ '/ 1 5 ;2 A c L /> * 

_ ^ ipc^.u'VAri.^.i 


133340 - 


-0 . 196273333+04 





r 


w 

w 




1 : 1 2 SB'S 


DATE 

SK : 

0.0000000 
0.0000000 
2.13 OOOOOE-02 

sc : 

0 . 0000000 
9 . 02999398-04 
ooooooo 


TIME 


b:io 


o. 
o 
y . 


ooooooo 
,0000000 
4349991 E-04 


0.0000000 
0.0000000 
1 . 46400008-03 


0.0000000 
0.0000000 
i .2i'jyyyyt-o3 


0. 
3 
0 . 


0 

so : 

0.0000000 
0.0000000 
4. 51 000008-04 

beta « o.ooooooo 
zer * -i.mttm-oj 
ak: 

2730000 . - 
-583.0310 
1 .82377998-02 
9 . 28 1 10028-0+ 

DM ««• -30505.23 

ROOTS : 

-0.831573911+02 
-0 . 8313739 IE +02 
-0.276121028+0 3 

- 0 . 7 7 6 1 7 1 0 2 E + 0 3 
-0.667192768+02 
-0. 3671927*8 + 02 

- 0 . 6 7 6 1 5 2 4 7 8 + 0 2 
476152 17E+02 
48321 190E+03 
43521 1908+03 
310525668+02 

-0.310525838+02 
-0.292751798+02 
-0.292751798+02 
-0.2038771 It +02 
-0.2048 77 11E+02 

- 0 . 2 1 5 7 5 8 4 B E + 0 2 
-0. 21 575848 E + 02 
-0.407595708+01 
-0. 407595708+01 
-0. 2572+9+8 t +02 
-0. 257 2 6948 E +02 
-0. 132890 + 9E + 02 
-0 . 132890+9E + 02 
-0. 10155215E+02 
-0 . 101552 15E + 02 
-0 . 87753722E + 02 
-0 . 877537228+02 

0.908925088+01 


OOOOOOO 
l 5000008-03 
OOOOOOO 


0.0000000 
2 .2789999E-03 
0.0000000 


0.0000000 
1 .87100008-03 
0.0000000 


0 . ooooooo 
0.0000000 
1.37600008-03 


0.0000000 

0.0000000 

1.09700008-03 


0.0000000 

0.0000000 

1.90000008-03 


2750000. 
-585.0510 
1 .82377998-02 
9 . 28 4 4 0028-06 


4050000. 
-814.8732 
-+. 38322838-02 
2.9026 4 00E-05 


4050000. 
-814.8732 
•6.3S32263E-02 
2. 90244008-05 


GAO« -0.2615030 


- 0 , 

- 0 . 

-0 

- 0 . 


0. 
-0 
0. 
-0 
0. 
-0 
0. 
-0 
0. 
-0 
0. 
-0 
0 . 


447321698+04 
41732169E+04 
409569348+04 
40 958934 E +04 
417151298+04 
41 715129E+0 | 
410415478+04 
410 41547 E +04 
341141348+04 
341141318+04 
334913438+04 
334915138104 
304590 1B8+04 


ORIGINAL 1 r i ^ 

OF POOR QO ™ 


-0 .30159018E+04 
0 . 2 y 3 9 6 i 6 9 8 + 0 4 
-0.293961898+04 
0. 27264 8 138+04 
-0 . 2 724 48 1 5E + 04 
0 . +877+0228+03 
-0 .687760228+03 
0.326969618+03 
-0 .52686981E+03 
0.275660778+03 
-0. 2 7 566077 E +03 
0 . 2 2 2 2 B B 9 3 8 + 0 4 
-0 . 2222B393E+04 
0 . 162B13448+04 
-0. 162315448+04 
0.206826738+04 


0 .90892503E+01 
-0. 2064 74 408+02 
-0 . 20817 4 408+02 
-0.974061878+01 
-0.974061378+01 


-0 .206326758+04 
0 . 132719348 + 04 
-0.182719348+04 
0. 193020028+04 


-0 


0 . 1 9502002E +0 1 


Bll 


J 


DA76 2J12J84 

SK : 

0 . 0000000 
0.0000000 
2.11000006-02 
sc : 

0.0000000 
?.o:»yyy«?K-oi 
0.0000000 
sc : 

0.0000000 

0.0000000 

4.51000006-04 


2730000. • 
-385,0310 
1 .82377996-02 
9.28110026-06 
0« "■ -30303.23 

f:outs : 

-0.J 20630916+03 
-0.120650916+03 
-0.638004896+02 
-0.438001396+02 
-0.417761398+02 
-0.417761596+02 
-0 . 239113038+03 
-0 . 239 1 130 38+03 
-0.483325738+03 
-0 . 1853257 5 8+03 
-0 . 390085868+07 
-0.390085868+02 
-0.29926 £588+ 02 
-0.299288386+02 


TIME 6 i 28 

0 . 0000000 
0.0000000 
9 . 43199918-04 

0.0000000 
3 . 13000006-03 
0.0000000 

0 . 0000000 
0.0000000 


2750000. 

-585.0310 

1.82377998-02 


0 . 0000000 
0.0000000 
1 .46400008-03 

0.0000000 

2.27899996-03 

0.0000000 

0.0000000 

0.0000000 

1.09700006-03 


4050000. 

-814.8732 

-6.38322836-02 

2.90264006-05 


9.28440026-06 

GAO" 0.2615050 

0.464278466+01 
-0.164278468+04 
0.418668756+04 
-0.418668756+04 
0.411020716+04 
-0.411020718+04 
0.391749646+04 
-0.391719646+04 
0.341111496+04 
-0.311111498+04 
0. •333484196 + 04 
-0.333484 196 + 04 
0.304067016+04 
-0.304067016+04 


1.37600006-02 
BETA « 0.0000000 

ZER <<• 4 /99999998-03 

AKt 


-0.233013136+02 

-0.235013136+02 

-0.22787406E+02 

-0.227874066+02 

-0.380171976+01 

-0.380171976+01 

-0.256521566+02 

-0.256521366+02 

-0.132820316+02 

-0.132820316+02 

-0.117914236+02 

-0.117914236+02 

-0.533217166+02 

-0.353217166+02 

-0.503762826+02 

-0.303762826+02 

0.252824626+01 

0.232824626+01 

-0.297373216+01 

-0.297373216+01 


0.291921336+04 
-0.291921336+04 
0.272659446+04 
-0 . 272659 4 4 F +04 
0.688235356+03 
-0.688235356+03 
0.326931526+03 
-0.524931326+03 
0.275664546+03 
-0.275664548+03 
0.221763706+04 
-0.221763708+04 
0.173797356+04 
-0.17379/336+04 
0. 184336726+04 
-0.181536726+04 
0. 193970196 + 04 
-0 . 1959/0196+04 
0. 192309S7E + 04 
-0 . 192309376+04 


OklGWAL PA. k ‘ .'J 
OF POOR QUALITY 

0.0000000 

0.0000000 

1.21399996-03 

0.0000000 

1.89100006-03 

0.0000000 

0.0000000 

0.0000000 

1.90000006-03 


4050000. 

-814.8732 

-6.38322836-02 

2.90264006-05 


812 







DATE 

sk : 
0. 
0. 


i : 12:34 


TIME 3 1 S 37 


0.0000000 0.00 
0.0000000 0.00 
2.11000008-02 9.43 

sc : 

0.0000000 0.00 
9 . 02WV9?8-0 4 n.i'j 
0.0000000 0.00 
so : 

0.0000000 0.00 

0.0000000 0.0'J 

4.53 000008-04 1.57 

SETA * 0.3000000 

ZER «• 4 . 99999998-03 

AK : 

2730000.- 275 

*999*0540 -58L 

1 .82377998-02 1.92 

9.29440028-06 ?*2« 

DM '* -30505. 2 3 

ROUTS : 

-0. 33 139733 8 +02 0.4 


0.0000000 

0.0000000 

y.434yyyit-04 

0.0000000 
3 . 1 300000+1-03 
0.0000000 

0.0000000 
0.0000000 
1 .37600008-03 


ORIGfNAL 
OF POOR 


0.0000000 

0.0000000 

1.46400008-03 

0.0000000 

2.278vyyyt-03 

0.0000000 

0.0000000 
0.0000000 
1 .09 700008-03 


PAGE (8 
QUALITY 

0.0000000 
0.0000000 
i .2i3yyyyt-03 

0 . 0000000 
1 . 89100008-03 
0.0000000 

0.0000000 
0.0000000 
V, 433001 Ok" -04 


2750000. 

-595*0540 


• 4050000. 
-814.8732 


1 .823779VE-02 -6 . 38322838-02 
9.28440028-06 2 . 9026 4008-03 

BAO» -0.2615050 


4050000. 
-814.8732 
-6.383 22 83E-02 
2.90264008-05 


447320828+04 
-0 .831397338+02 -0 . 4 47320828 + 04 
-0. 2 7 614064 F +03 0.409585138+04 
-0.276140648+03 -0.409385138+04 
-0.677898848+02 0.410413448+04 
-0 .677898348 + 02 -0.410415448 + 04 
-0 . 867081928 + 02 0.41 715073E+04 
-0.967081928+02 -0 . 417130758+04 
-0.483136068+03 0 . 34 1 1 34 668+04 
-0 . 433136068 + 03 -0 . 3 41 1 34668 + 04 
-0.341073518 + 02 0 . 334921 208+04 
-0.341073318+02 -0 .334921208+04 
-0.292686758+02 0 . 304 389 1 38+04 
-0. 292686 738 + 02 -0.30 4589138 + 04 
-0.205222428 + 02 0 . 293965208" + 04 
-0.20322242 8 +02 -0.29396 5 208 + 04 
-0.213293408+02 0.272631598+04 
-0.215293408+02 -0.272651398+04 


-0.409176498+01 

-0.409176498+01 

-0.237372718+02 


0.687784698+03 

-0.687784698+03 

0.526878708+03 

-0.526878708+03 


-0 .25/372718+02 -0 . 526878708 + 03 
-0. 132886588+02 0.273638348+03 
-0 .132886388+02 -0 . 273638348 + 03 
-0.897806658+02 0.162342878+04 
-0.897806638+02 -0 . 162342878+04 
-0.100203328+02 0 . 222282588+04 


-0 .100203328+02 -0.222232338 + 04 
-0.2097841SE+02 0.182730948+04 

-0.209736158+02 -0.132730948+04 
0.110285018+02 0.206928318+04 

0 .110235018+02 -0.206323318 + 04 


-0 .968613678+01 
-0.969613678+01 


0. 193022268 + 04 
-0.193022268+04 


v=) 


B 14 
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D.4TE 2:iO:E4 

TIME 15 5 0 

OF POOR QUALITY 

sk: 




0.0000000 

0.0000000 

0.0000000 

0.0000000 

0.0000000 

0.0000000 

0.0000000 

0.0000000 

2.110 0 0 0 0 E •• 0 2 

C » 

9. 434999 IE-04 

1 . 4£ 49009E-0S 

l . 2159999E-03 

0.0000000 

0.0000000 

0.0090000 

0 * OOCOOO 0 

9 .0299927E-04 

3 . 1300000S-03 

~ n-»-3oooaiT-A-r 

•>>*./ W / ( / / w */w 

1 . 39 1 0000E-C3 

0.0000000 

0 . 00 00000 

0 # 0 0 0 0 0 0 0 

o . ocoocoo 

ca: 




0.0000000 

0.0000000 

O.OOCOOCO 

0.0000000 

0. 0000000 

0.0000000 

0.0000000 

0.0000000 

4 . 5 1 OOOOOE-04 

1 .57600C9E-03 

1 .0970000E-03 

9.4520010E-04 

9ETA = 0.500000 




ZER = 4 . 9999999E-02 



ak: 




2750000 . 

2750000. 

4050000. 

4050000. 

-533.054'.'' 

-525.0540 

-314.372? 

-314.9732 

1 • 223779? S> 02 

1 .E22779 7E- 02 - 

£ .38322838-02 

-5 . 3P222S2E-02 

9.2344002E-06 

9*28440028-01 

2*90264008-05 

2.9C26400E-C3 

DM = -20305,23 

ChC= 0*261 

5050 


roots: 




-0.12060695E+03 

0. 46427ES5E+04 



-0 . 120£0£95E+02 

-0.46427233E+04 



-0.637S7359E+02 

0.41S££3?7E+04 



-0 .62727359E1-02 

-0#413563?7E*01 



-0 . 418339798+02 

0.41102093E+04 



-0 • £ 1 2337702 f 02 

-0.411020932+04 



-0 . 21 941 1 9GE+03 

0*391750648+04 



-0 .2594519CE +03 

-0*27175 034S+C4 



. 4S2 1 752 ££t 03 

0#34110472£+01 



-0.42S17526E+02 

-0.341104728+0 1 



-0 . 3915S422E+02 

0 . 333482838+0 ' 



— A lOlHAmciA') 

v * w « i k. W 7 . w W. 1 w *_ 

-0*333432838+04 



-0*300837808+02 

0 . 30406552E+0 » 


# 

-0.30023730E+02 

-0 .304065E9E+04 



-0*234078038+02 

0.29192920E+04 



-0.234072 9 3E+02 

-0.2919293CE+04 



-0.2276S338E+02 

0.27265741E+04 



-0.22732333Et-02 

-0 . 27265741E+04 


# 

-0 . 3 30 7 1221 £+01 

0, ££32524 0E+03 



-0*380712218+91 

-0 . 63825240E+03 



-0 • 2545 99S4E+0? 

0 . 52692522E+02 



-0 *250599548+02 





-0. 132S2979E+02 
-0 . 1 32929 79E+02 
-0. 11719 452E+02 
-0 . 1175 o 452E+02 
-0.564622 77E+02 
-0 .56462277E+02 

_ A . «■» 0?AAO-7CJ.A-> 
-n . 1 A7/.AO-?Ci,V) 

A T0JO7O1 

«-• ■/ W . a W ■» w - » 


0.27566517E+02 
-0 .27566517E+03 
0.22175953E+04 
•C . 22175953E+04 
0. 173920S9E + C4 
■0 . 17392059E+04 
0 . l£444£05Z+0* 
■ 0 . 1344 4 £051+0 •* 

0 . 1 « ^ 51 £ £ i — A A. J 


rv tiooioinriAi 

W • WrfW W W A W la > t/ * 

^ 4 C « IA1 

• *; * 2 ! 5 i A 1 2 ,r L 0 1 


1 ° 5 ? 1 CCwE r 

1 ° 2 3 o 1 7 ° E *f* 0 • 
- 1 ? 2 2 £ \ 7 ? 3 ^ 0 •* 



ORIGINAL PAGiT »3 
OF POOR QUALITY 


Channel 24 


Channel 11 


Recorder 2 
2 vol ts/1 ine 
.5 mm/sec 


Recorder 1 
2 vol ts/1 ine 
.2 mm/ sec 


CDCZ 6-27-83 
TWOZETR = .3 





« 


i : j 1 ; 8 4 


DATE 

sk : 

0.0000000 
0.0000000 
2. 11000008-02 
sc : 

0 . 0000000 
9 . 02999398-01 
0 . 0000000 
su : 

0.0000000 
0.0000000 
4.1=13 000008.-04 
BETA '< 0.0000000 

ZER «• 0.11500000 


TIME 1 4 5 33 


ORIGINAL ? 

or foor Q 


\‘;V. (3 

jaLITV 


0.0000000 
0 . 0000000 
9.43499918-04 


0.0000000 
0.0000000 
1 .46400008"-0i5 


0.0000000 
0.0000000 
1 .2159999E-03 


0.0000000 
3 . 10000008-03 
0.0000000 


0.0000000 
2 . 27H9999E-03 
0 . 0000000 


0.0000000 
1 .8910000E-03 
0.0000000 


0.0000000 
0.0000000 
1 .5760000E-03 


0.0000000 
0.0000000 
1 .0970000E-03 


0.0000000 
0.0000000 
9 » 4530010E-04 


ak: 

271=50000. 
-380.0340 
1 .82377998-02 
9.234 10028-06 
DM "• -30S03.23 

ROOTS : 

-0. 1786161 8E+ 03 
-0. 178616138+03 
-0. 63234 4608" + 03 
-0 . 652344608+03 
-0 ♦ 867326328" + 02 
-0.867326528+02 
-0. J 1363S63E+03 
-0.113635638+03 
-0.488000532E + 03 
-0.188000528+03 
-0.90215115978 + 02 
-0.9025t597E+02 

- 0 . 3 8 4 8 1 4 3 1 E + 0 2 
-0.38431 4318+02 
-0 ♦ ’367208 12E + 02 
-0.367208128+02 
-0.8911 2V4SE+02 
-0.891 129408+02 
-0.382552708+01 
-0 . 332552708+01 

- 0.2=31501 8 02E +02 
-0. 20001 802E+02 
-0 . 5 32639078 + 02 
-0. 132639078+02 
-0 . 3 3229742E+02 
-0.132297428+02 
-0.230674248+03 
-0.230674248+03 

0 . 6389077*3E + 0 1 


2700000. 
-080.0040 
1 . 8 2 3 7 7 9 9 E - 0 2 
9.23440028-06 


4000000. 
-814.0732 
- 6 . 3 8 3 2 2 8 3 E - 0 2 - 
2.90264008-00 


4000000. 
-314.3732 
6 . 3832283E-02 
2.90264008-05 


DAO--- -0.2610000 


0. 44497290E+04 
-0.444972908+04 
0 ♦ 3 9 8 0 1 0 2 3 E + 0 4 
-0.398010238+04 
0. 4 17208038 +04 
-0.417208038+04 
0. 40802 646E + 04 
-0.408026468+04 
0 . 3 4 1 0 8 4 27 E + 0 4 
-0.341034278+04 
0. 336321 00E + 04 
-0.336321008+04 
0.30444089E+04 
-0.304440398+04 
0.293888778+04 
-0.293838778+04 
0.27633658E+04 
-0.276336088+04 
0 . 6 8 7 7 0 2 8 5 E + 0 3 
-0.687702308+03 
0. 32680899 E +03 
-0.026808998+03 
0. 27067031 E+03 
-0.270670318+03 
0.2 2230390 E+04 
-0.222300908+04 
0.1 71080938" + 04 
-0.171030938+04 
0 .20701221 8 + 04 


0 . 638907708 + 0 1 -0 . 20 70 12218 + 04 
-0.38167225E + 02 0 . 1 80739 1287 + 04 

-0.331672208+02 -0.130/39128+04 
-0. 110614 49E + 02 0 . 192921 368" + 04 

-0.11061 4 198+02 -0.192921368+04 


4 


B17 


* 4-1 



4 

I 


DATE 2:ii:81 
SK : 

C. 0000000 
0.0000000 
2 . 1 J OOOOOE-02 

sc : 


ORIGINAL PAC.L .3 
TIME 11:12 OF POOR QUALITY 


0 . 0000000 
0.0000000 
y . 4349991 e-oi 


0.0000000 

0.0000000 

1 . 4640000E-05 


0.0000000 
0 . 0000000 
i . 21399993-03 


0.0000000 
9 . 02999398-01 
0.0000000 
su : 

0 . 0000000 
0.0000000 
4«9SOOOO0E-O4 
BETA * 


0.0000000 
3 . 1500000E-03 
0.0000000 

0.0000000 

0.0000000 


0.0000000 

2.278vyyyt-03 

0.0000000 

0.0000000 
0.0000000 
1 . 0970000 E -03 


1 .37400003-03 
0.0000000 


ak: 

2750000 . 

-585 . 0 Li 1C 
1 .8237799E-02 
9.2844OO2E-06 
DM <• -305013.23 

ROOTS : 

-0.22141724E+03 
-0.221 11 724E+03 
-0.48195993E+02 
-0.48193993E+02 


2750000, 1050000. 

-383.0340 -S1 1 .8732 

1 . 8237799 E-02 -6.3832283E-02 
f. 284400 2E -04 2 . 9 0 2 6 1 0 0 E - 0 8 

800" 0.2615050 


0. 44241491E + 04 
-0.44241491E+04 
0. 41440 7 81 E +04 
-0.41440781E+04 


-0.3 1 707 1 20b +03 
-0 » 11707120E + 03 
-0. 42304 782E+03 
-0,423 0178 2E+03 


o.iioyosy?t:+04 
-0 . 410908892 + 04 
0.374439088+04 
-0 .37443908E+04 


-0.48920330E+03 
-0.43920350E+03 
-0.1 1577530E+03 
-0. 1 1377530E + 03 
-0.43497914 E +02 
-0. 4349791 4E +02 
-0.990341 29E+02 
-0.9903 4 129E+02 
-0. 7174081 1E+02 
-0. 717 4 0811 E + 02 
-0 . 2363 1 734 E + 02 
-0.2343 1731 E +02 
-0.37217203E+01 
-0 .37217203E+01 
-0 • 1 32879 8 9 E+ 02 
-0 » 1 3287989E+02 
-0. 14912758E+02 
-0.14912738E+02 


0 . 3 4 1 9 4 5 4 5 E + 0 4 
-0.31194545E+04 
0, 3 3 43821 7 E +04 
-0.33458217E+04 
0. 30371 008K+04 
- 0 . 3 0 3 7 1 0 0 8 E + 0 '» 
0 . 293 4 54 80 E+04 
-0 .2934348CE+04 
0. 27842421 E+04 
-0.27812421 E+04 
0.52491252E+03 
-0.S2491253E+03 
0.4B815943E+03 
-0.4S815913E+03 
0 , 275474 22E+03 
-0 .27547422E+03 
0.221 78381 E+04 
-0 .22 173381 E+04 


0.0000000 
1 .3910000E-03 
0.0000000 

0.0000000 
0.0000000 
9 . 45300 1 0E-04 


4050000. 
-814.8732 
-4. 3 83228 3 E! - 02 
2 , 9 0 2 4 4 0 0 E - 0 5 


-0.31V14394E+02 0 . 177451 80E+04 
-0.31914394E+02 -0 . 17713130E+04 
-0. 19108209E+03 0 . 1 884 1 393E +04 
-0 » 19 10820 9 E +03 -0 . 18841898E +04 


-P-t-l -liP 2 2 4 4 4 E + 0 2 
0.15022 14 4E+02 
-0.4325024 If +01 
-0 . 13250211E + 01 


0 .19575421 E+04 
-0 .195754218+04 
0. 1925 87 OOf +01 
-0, 192337 COE+ 04 


I 






\ 

i 

> 

1 

1 


* 

1 



bis 


J 
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APPENDIX C FOURIER FILTER FORTRAN CODE 


DP 1 *BFFIL . FOR JANUARY 16,1934 1^:00 

FOURIER TRANSFORM FILTER FOP HYBRID SIMULATION SINUlkTIvn 


H T U T VARIABLES 


FILFRE 

Filter f r e q u e r. c y 

FPEINT 

Frequency i ntirva 1 

DT 

S i iTi u I a t i c n t i m i stc p 

TFINAL 

Stop t im-: 

NPTS 

Number of si mu 1 1 1 i on po i nt / 

RNFIN 

N u m b s r o f F o u r i € r i n t « r val • 

ISPT 

Sample c v e r y 1 3 PI point:* 

NP7PLT 

Plot every HP7PLT point/. 

JN • 

Joint number 



INTERMEDIATE 

variables : 


NN 

Number of point/ in Four id 

C 

K'PRIME 

FILFRE hj s • / e n u n b e r 

A V 

K 1 

First wavenumber of freque 


KF 

L a j t w a v e n u m b e r of f r e q u e n 

t' 

OMEGAK< 9 ) 

Fundamental frequency for 

c 

COS 1 K< 9 > 

C o rind o f C M EGA K 


SIN1 K< 5 > 

Sine of GMEGAK 

1 . 

TO 

Initial tin * of z i m u 1 a t i o n 

c 

I PT 

Humber of data point/ used 

f 

IFTFLT 

Plot point counter 


OFT 

F 0 counter 

c 

IPLFT 

N u n b e r o f plot p o i n t s 


OUTPUT VARIAB 

LES : 


T 

Time 


AK< 9 ) 
BK< 9 ) 

F. M A G < 9 > 


F o u r i e r c o eft" i c i e n t s 
Fourier coefficients 
M a 0 n i t u - J c o f A K and B t . 


t > I M E M 3 1 0 N A K \ 9 , B K ( 9 > , 0 M EGA K < ? > , R M A G < 9 *> , F < 5 2 0 C > 
D I MENS I ON COS 1 K< 9 > , S I H 1 K< 9 > 

CHAR ACTER+2 1 I F M T 

ccept par abaters from term i n -a 1 

WRITE'-' 7,*> 'Input frequency of interest ( r ad/se< 
READ< 3 , * ) FILFRE 

WRITE<7, + > 'Input frequence interval' 

READ<5,+) FREINT 

WRITE'?,*) 'Input simulation time interval' 

RSAD< 3 , * ) DT 

WRITE-: 7, *> Input final time' 

READ-! 3 , * > Tr INAL 

w " I T E' 7 , * > 'Input nuffiwsr of points to consider 
a E n D > 5 . * » h - T 

WRITE 1 ' 7 , *> 'Input nun ter o** Fourier intervals tc 
READ'' 3, •*■ ) RN’F I M 


Cl 


7 11:23 1334 


bffil.f Page 


WRITE<7,*> 'Sample every n points; input n ' 
READ<5,*> I3PT 

U?ITE<7,*> 'Input points/plot ' 


OF 1 ? L PAC * » 

p OOR QUALITY 


REhD< 5 , * ) NR TFLT 

WRITEC 7, + > 'Input Joint number -3 1 sp 1 segments ' 
JsEAD<5,*> JH 


$.» t up variable format statement to read every I 5PT points 

I F< I5PT.EQ. 1 ) THEM 
WRITES I FMT / 6 0 0 0) JN*1 1 


ELSE 

WR I TE< IFMT,70Q0> I5PT-1 , JN + 1 I 
EMC- IF 

0 F0RMAT< '< ' , I 3 > ' < / > , ' , 1 3, 'X , E 1 1 . 4 > ' > 
0 FORMATS '< ' , I 3 , 'X , E 11 . 4 ) ' ) 


Compute filter time interval 
>T * D T + 1 S P T 


Opert file where initial time and DtLY data are stored 

OPEN < UN I T= 1 0, FILE= 'H YE 1 .DAT', STATUS* 'OLD ' , FORM* ' FORMATTED ' > 

Open files to send plotting data 

OPEN < UNIT *13, FILE* ' 6 F F I L I .OAT ' , F CRM* 'FORMATTED ' > 

OPEN < UNIT* 15, FILE* 'BF FILS .DAT' . FORM* FORMATTED ' > 

OPEN < UNIT = 1 7, FILE* 'EFFIL5 . DAT ' .FORM* 'FORMATTED ' ) 

C a leu! at e number of pc ints necess ary f or each Four ier per i od and 
frequency interval 

NN*IF IX< .5+2. *3 . 1 4 1 53265/C DT*FRE IHT ) > 

FREINT*2 . *3 . 14159265/< NN + DT > 

I TEMP* I F IX< RHF I N + HN > 

IF< ITcMP.LT.NPTS) NPTS* I TEMP . 

IF< NN . LE . 520 0 ) GOTO 7 0 

WRITE*' 7, + ) 'Number of points in Fourier interval is greater 

URITEC7,*) 'dimension of F NN* ' , NN 

STOP 

Determine K values for region of interest 

Dimensions of AK, BK, QMEGAK, RMAG, COS 1 K, and SIN IK must be changed 
number of K's is increased 

KPRIHE*IFIX< . 5+FILFRE/FREINT > 

K 1 *KPR IME-4 
NK*9 

IFCK1 . GT . 0 )G0T0 SO 
K1 *1 

NK-*K.--K1 + 1 
KF=KPR IME+4 
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ORIGINAL PAGE 19 

OF POOR QUALITY 


Zero arrays: calculate constants for each frequency 
J a 1 

DO 10 K«K1 , KF 

OHECAK< J J-FLOATC K >+f REINT 
COSt K< J )«C0S< OMEGA! < J >*D7 > 

SIH1K< J >»SIN< GM£GAK< J > * D T > 

AK< J >*0 . 

BK< J>«0. 

•)- J+ 1 

Initialize time, point counter, plot counter, and F0 pointer 

F.EAD< 1 0, 1 000 > T 
T0 = T 
I FT = 1 

IPTF'LT =MPTFLT 
IFLF'T = 0 
Jr T = 1 

M a i n 1 o o p 
CONTINUE 

Read DELY = FCn+1 ) 

READC 1 0, IFMT > FNP1 


Adjust counter to proper F 0, get F0 and replace with F< r* J > 

IFf. JRT , GT . < MN+ 1 > ) JF 7= JPT-MN- 1 
F0 = F< JPT > 

F< JPT )*FNP 1 

Calculate coefficients 
J=1 

DO 3 0 K * K 1 , K F 

TEMPI *AK< J >*S I N 1 K< J ) 

AK< J )*AK< J >*COS1K< J > + BK< J >*SINi K< .J >-FG*C0S 1 K< J 3+FNF1 
BK< J >*F 0+SIN 1 K< J 3-TEMF 1 +BK< J>*COS1K< J> 

RMAGC J >=SQRT< AKC J >**2+BK< J >**2 ) 

J= J+1 

Test for output conditions, write time and coefficients to files, 
increase plot counter 

Ir< IPTF'LT. NE.NPTPLT) GOTO 90 
WR I TE< 13,3000 ) T , < AK<. J > , J« 1 , nk > 

WP I TE< 1 5, 300 0 > T , < BKC J >, J=1 , nk > 

UR I TS< 1 7, 300 0 > T , < RMrtGC J ), J* ! , nk ) 

IPTPLT*0 
IPLPT-IPLPT+1 
I:-7FLT = IPTrLT + 1 

Check termination conditions 
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OR'TNV. FAGr I> 

Of- FOUR QUALliV 


IF< <. T . CE . TF IHAL > . OR . < IPT . GE . HPT S ) > GOTO 50 
Increment t if-, and point counters 


T = T T 
IPT«IPT*1 
JPT*JPT+1 
GOTO 100 
CONTINUE 


Close output files 


CL03EC UH I T* 1 0 ) 
CLOSEC UNIT* 1 "3 ) 
CLOSEC UN I T* 1 5 ) 
CLOSEC UH I T = 1 7 ,* 



Op -in files for plotting data 

OPEN < UN IT=14, FILE*' EFFIL2.DAT' > 
OPEN C UN I T* 1 6 , F ILE* 'BFF I L4 . DAT ' ) 
OPEN < UN I T* 1 8 , F I LE= 'BFFIL6.DAT' > 


Sind number of points and variables to plotting files 


WRITEC 14,3000 > 10,IFLPT 
WRITEC 1 € ,5000 > 1 0 , IPLRT 
WRITE*'. 13,3000) 10,IFLFT 

Write titles for plots 


WRITE-'. 14,6000 ) 

'TIME' 

WPITEC 16,6000) 

'TIME ' 

WRITE 1 ' 13,6000 > 

'TIME ' 

WR I TEC 14,6000 ) 

' AKC 1 > ' 

WRITE*' 1 4,6000 ) 

'AKC 2 )' 

WR I TE< 14,6000) 

'AKC 7 * 

WRITEC 1 4,6000 > 

'ARC 4 ) ' 

WRITEC 14/ 600 0 > 

'AK* 5 ' 

WRITEC 14,6000) 

'ARC 6 > ' 

WRITE*' 1 4,6 000 > 

'ARC 7 ) ' 

WRITEC 14 6000) 

'ARC 3 > ' 

WRITE*' 1 4 ,6000 > 

'ARC 9 ) ' 

WRITEC 16,6000 ) 

'BK< 1 >' 

UR I TEC 1 6, 6000 > 

'BKC2 >' 

WRITEC 16,6000 ) 

' BKC 3 ) ' 

W R I T E C 16,60 0 0 ) 

' BKC 4 ) ' 

WRITEC 16,6000 ) 

'BKC 5 > ' 

WRITE' 1 6,6000 > 

'BKC 6 >' 

UR X TEC 1 6 • 6000 > 

'BKC 7 > ' 

WRITE-: 16, i 0 3 ) :• 

* z V , • ‘ * 

Ui. \ w . 

WRITEC 16 , 6 0 0 0 ) 

'BKC 9)' 


-V 
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UR I TEC 13,6000) *FMAG< 1 > ' 

WP i TE< 13,6000) 'FflAG'2)' 

WRIT£< 13,6000) 'Rt1AC< 3 ) ' 

WRITE< 13,6000 'RHAG<4>' 

WR I TE< 1 3 , 6 0 0 0 ) ' F MAG*' 5 > ' 

* URITE< 13,6000) 'RMAC'6>' 

wr< ! TE< 1 3 , 6 0 0 0 ) 'RMhG 1 ' 7 > ' 

WRITE-: 13,6000 ) 'FMAG' 3 > ' 

WR I TE< 13,6000) 'PMAG< 9 ) ' 

C 1 © s -s p 1 o t *• i n f i 1*5 

CLOSE* UN IT* 14 ) 

CLOSE*: UNIT* 16 > 

CL0SE< UNI T* 1 3 ) 

C Writ* piriinitirf to fili 

OPEN*: UNIT* 19, FILE* 'EFFILF , DAT ' ) 

UP. I TE< 1 9 , * ) ' NN : ' , NN 

WRITE' 19,*) ' D T : ' , D T 

WF ITE< 19,*) ' FILTER FREQUENCY j ' , F I LFRE 

WR I TE< 1 9 , * > ' FREQUENCY INTERVAL : ' , FFE I NT 

WRITE<19,*> ' DATA GOOD AT T=',T0 + DT*FL0AT< NN ) 

WRITE* - .19, * ) ' JOINT NUMBER : ' , JN 

WRITE*: 19, * ) 

WR ITE*: 1 9, * ) ' FCN K Cr'EGHr 

WRITE-' 19,*) 

J=t 

DO 4 0 K = K 1 , K r 

W= ITE' 19,2000 > J , K , 0MEGAK( J > 

0 J = .J- 1 

CLOSE' UN IT* 19) 

C 

Fermat 5 tat*=ni*r.t5 

1 000 FORMAT* El 1 .4 > 

7 000 FORMAT*: 1 0E1 1 .4 ) 

000 FORMAT*: 21 5) 

2000 FOPMAT<SX, 12, 4X, I4,SX,F1 1.3) 

60 0 0 FORMAT-: AS) 

END 
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HYBRID SIMULATION VARIABLES 

The need for better documentation became apparent when we began to study the 

A 

hybrid simulation code in depth. An outdated list of variables was expanded 
(Appendix D) , a list of comparable stability model variables compiled (Appendix 
E), and a comparison list made (Appendix F). These lists are not meant to be 
totally inclusive in themselves but provide brief descriptions of the important 
variables and some minor ones. 

Where there is a common variable at each of the seven interaction joints it 
is listed as VAR(P-S3), where VAR represents the part of the variable in common. 

For example, the variable listed DEL(P-S3)(Y,Z) includes DELPY, DELPZ, DELTY, 

DELTZ, ... for P, T, B, A, SI, S2, and S3 interaction joints (see Figure 4, 

Section 2.3). 
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HYBRID SIMULATION VARIABLES 


Description 


Simulation Variable 


Case accelerometer outputs 

Accelerometer modal coefficient 

Balance piston damping curvefit coefficients 

Balance piston spring curvefit coefficients 

Bearing or bearing + carrier (dependent on 
IBSEPBC) stiffness coefficients between deadbands 

Curve fit coefficients for AK(P,T)(Y,Z) 

Bearing stiffness coefficients inside deadbands 

Bearing stiffness coefficients outside deadbands 

Principal inertia about (X,Y,Z) axis for Jth 
rotor joint 

Mass of Jth rotor joint 

Type of pump for output only 

Scale factor for Alford cross-coupling 
coefficient 

Stiffness of backup structure at bearings 

Buffers to hold channel outputs 

Damping coefficients 

Squeeze film damper coefficient 

Squeeze film damper coefficient 

Bearing damping coefficients between deadbands 

Bearing damping coefficients outside deadbands 

Bearing damping coefficients inside deadbands 

Multiplication matrix for deflections at Jth 
accelerometers 

GP(P,T) - G(P,T) constant 


A( 1-7) 

AC( 1-7) 

ACB(0, 1 ,2) 

AKB (0,1,2) 
AK(P,T)(Y,Z) 

AK(P,T)( Y.Z) (0-4) 
AK(P,T)DB( Y,Z) 
AKP(P,T)(Y,Z) 
AIR(1,2,3)(J) 

AMR(J) 

APUMP 

BETA 

BK(P,T)(Y,Z) 
BUFFR(A.B) 
C(B,SA,S1,S2,S3) 
C( 1,2,3) 

C( 1,2,3) (P,T) 
C(P,T)(Y,Z) 
CP(P,T)(Y,Z) 
C(P,T)DB(Y,Z) 

CAC( 1-7) (J) 

CL(P,T) 
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HYBRID SIMULATION VARIABLES, con’t 


Description 


Simulation Variable 


Center of mass location 
Cos(PHIX) 

Output to channel I 

Maximum value during data sampling for channel I 

Minimum value during data sampling for channel I 

Average DA1MIN over data gathering Interval 
Average DA1MAX over data gathering interval 
Entire selection of output values 
Accelerometer position relative to joint J 
Derivative of DPP 
Absolute magnitude of DEL 

Relative displacement from centerline between 
rotor and case 

DEL in X direction at balance pistons 
Derivative of DEL 
Derivative of DELBX 
Derivative of DELP 

DEL between inner and outer deadbands 

DEL with outer deadband removed 

Maximum dimension considered to be zero 

Determines PHIDDX for ramp runs 

Integration time step 

DT/2 

DT**2/4 

Program default DT 

Input length for calculating C(1,2,3)(P,T) 


CM 

CX 

DAI ( I ) 

DAIMAX(I) 

DAIMIN(I) 

DAIMNAVG(I) 

DAIMXAVG(I) 

DAC 

DAC( 1-7) ( J) 
DDPP(P,T)( Y,Z) 
DEL(P,T,A,S1,S2,S3). 
DEL(P-S3)(Y,Z) 

DELBX 

DELD(P-S3)(Y,Z) 

DELDBX 

DELDP(P,T)( Y,Z) 
DELP(P,T)(Y,Z) 
DPP(P,T)( Y,Z) 

DMIN 

DSPEED 

DT 

DT02 

DT04 

DT2M15 

EL(P,T) 
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HYBRID SIMULATION VARIABLES, con't 


Description 


Simulation Variable 


Time elapsed while gathering data ELAPTIME 

C3(P,T)**2 EL0D2(P,T) 

Pump and turbine max forces when cavitation is EMAX(P,T) 
used 

Normalized case bending modal displacement ETA,ETAD,ETADD(J) 

and derivatives 


ETA,ETAO,ETADD of T - DT 
Not used 

Rotor forces except imbalances 
Rotor axial force 
Force used on line 734 
Case mode normalized translation 
Rotor mode normalized translation 
Total forces applied to rotor 
Maximum stiction level at SI 
Bearing squeeze film forces 
Pump imbalance forces 
Rotor modal mass integrals 
Pump bearing clearances 
Turbine bearing clearances 
Factor for DSPEED 
Modal mass integral 

Logical variable indicating position of switch 
to reverse spin speed derivative 

Flag: if * 1, "stiffness constants for bearings 
and their respective carriers are to be consid- 
ered independently" ; if = 0 , dependent 


ETAB,ETADB,ETADDB( J) 

ETABB.ETADBB, 

ETADDBB(J) 

F ( P— S3) ( Y , Z ) 

FBX 

FCAP( P,T) ( Y,Z) 
FEEC(I.J.K) 
FEER(I.J.K) 
FR(P-S3)(Y,Z) 

FS1MAX 

FSF(P,T)(Y,Z) 

FU( P-S3) ( Y.Z) 
GAM(O.l) 

GP.GPP 

GT.GPT 

GSPEED 

1 ( 1 , 2 ) 

IBACK 

IBSEPBC 


Case connection indices 


IC(P-S3) 
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HYBRID SIMULATION VARIABLES, con't 


Description 


Simulation Variable 


Squeeze film damping flag ICAV 

(O=none,l=cavit ,2=non-cavit) 


Name and dimensions of channel assignments 

Switch position for updating files 

Provision for modifying channel assignments 

Switch position to exit program 

Flag to indicate inclusion of bending modes 

Switch position for gathering output to 
line printer 

Maximum number of passes before buffout 
IMAX x 10 

Index for output to tape DAC channel 
Rotor connection indices 

Logical variable indicating position of switch 
for run time input 

Flag to signal Seal 3 inclusion 

Switch number of ITAPE 

Switch position for tape update 

Sideload flag (-I=zero,l=recalculated,0=unchanged) 

Marder data word 

Joint locations for accelerometers 

Number of case joints 

Number of rotor joints 

Flag to indicate update data case 

Stiffness coefficients 

Bearing stiffness coefficients between deadbands 
Number of data minus one to be output to tape 
Bearing stiffness xoef ficients outside deadbands 


ICARD(I.J) 

ICARDS 

ICHNNL(K) 

I EX I T 
IFLEX 
IGATHER 

IMAX 

IMAX10 

IOUT(K) 

IR(P-S3) 

IREMOD 

IRS3 

ISLTAP 

ITAPE 

ISIDEL 

IWORD 

JAC( 1-7) 

JNTC 

JNTR 

JREMOD 

K(B ,A,S1 ,S2,S3) 

K(P,T)(Y,Z) 

KCHAN 

KP(P,T)(Y,Z) 
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HYBRID SIMULATION VARIABLES, con't 


Description 


Simulation Variable 


No of passes (integration steps) during data 
collection 

KGPASS 

Flag: 1 Invalid case; 2 HPFTP with PH I X>0 ; 
3 HPOTP with PHIX<0 

KPUMP 

Position of interaction joints along rotor 

L(P,T,A,B,S1,S2 

Position from center of mass along rotor in 
x direction of joint J 

LR(J) 

Total mass of rotor 

M 

Maximum no of int. steps while collecting data 

MAXGPASS 

Number of case modes 

MO DC 

Number of rotor modes 

MODR 

Max number of passes before buff-out occurs 

NMAX 

Number of passes 

NPASS 

Run number 

NUMRUN 

Case mode natural frequency of Kth mode 

OHMC(K) 

0HMC( K)**2 

OHMCSQ 

Rotor mode natural frequency 

OHMR 

0HMR**2 

OHMRSQ 

Rigid rotor rotation 

PHI(X.Y.Z) 

Case mode shape translations (before assignments) 

PHIC 

Case mode shape translations 

PHIC(P-S3)(Y,Z) 

Case mode shape translations 

PHIC( P,T,B) X 

Rigid rotor angular velocity 

PHID(X,Y,Z) 

PHID of T - DT 

PHID( X, Y,Z) B 

Rigid rotor angular acceleration 

PHIDD(X,Y,Z) 

PHIDD of T - DT 

PH1DD(X,Y,Z)B 

Absolute value of PHIDX 

PHIt'XA 
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HYBRID SIMULATION VARIABLES, con't 


Description 

Simulation Variable 

PHIDX**2 

PHIDXSQ 

Rotor mode shapes before assignments 

PHIR 

Rotor mode shapes 

PHIR(P-S3) 

•R/2 

PI2 

Squeeze film damping coefficients 

PSF(P.T) 

Rigid rotor modal rotation 

PSI(X.Y.Z) 

Local rigid case modal rotation 

PSIC( A,B) X 

Cross coupling coefficients 

Q(A,SA,S1,S2,S3) 

Alford force cross-coupling curve fit coefficients QA ( 0 - 4 ) 

X axis location of Jth joint 

R(J) 

Squeeze film damper parameter 

R(P,T) 

FS1/FS1MAX 

RS7 

Seal damping curvefit coefficients 

SC(A,1,2,3)(0,1,2) 

Scale of Kth recorder channel 

SCALE(K) 

Full scale of Kth recorder channel 

SFI(K) 

Seal spring curvefit coefficients 

SK(A,1,2,3)(0,1,2) 

Sideload coefficients for straight line 
7 segment curve fit at balance pistons 

SL(Y,Z) 

SL at seals 

SL(A,S1,S2,S3)(Y,Z) 

Sideload force factors (dependent on ISIDEL) 
at balance pistons 

SLK 

SLK at seals 


Sideload forces at balance piston 

SL(Y,Z)0 

SLO at seals 

SL(A,S1 ,S2,S3) ( Y,Z)0 

Switch number for runtime input 

SLRMOD 

Intermediate coefficients in computation 
of sideload forces 

SL(Y,Z)T 

SL(A,S1,S2,S3) ( Y,Z)T 
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HYBRID SIMULATION VARIABLES, con't 


Description 


Simulation Variable 


Rotor speed time line factor 

Seal cross-coupling curvefit coefficients 

S I N ( PH IX) 

Sum of torques acting on rotor 
Intermediate sideload factor 
Final time 

Direction of sideload forces 
Direction of sideload forces 
Sideload phase angle 
Rotor speed time line factor 
Rotor speed at Kth mode 
Tic mark spacing 

2*DT; increment for speed reversal (IBACK) 

Damping coefficient of case mode K 

Damping coefficient of rotor mode * 

(GP(P.T) - G(P,T))/DEL(P,T) 

Fluid viscosity for squeeze film damping 
coefficients 


SPEED(K) 

SQ(A,1,2,3)(0,1,2) 

SX 

TAUM 

TEMP 

TF 

THSL 

THSL(A,S1,S2,S3) 

THSLOS 

TIME(K) 

TPHIDX(K) 

TPULSE 

TWODT 

TWOZETC(K) 

TWOZETR 

U(P.T) 

UM(P,T) 


Rotor imbalance moments 

(DEL(P.T) - GP(P,T))/DEL(P,T) 

Whirl orbital angular velocities 

Rotor center of mass X displacement and deri vativesX.XD.XDD 

XD.XDD of T-DT XDB,XDDB 

Nominal units/line on recorder XF1 

Rotor coulomb friction model denominator minimum XIDO 
velocity 


UM(P,T,B,A,1,2,3)(Y,Z) 

UP(P.T) 

VW(T,P) 


ABS(XTEMP)/ XIDO 


XIDMAG 
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APPENDIX E 


STABILITY MODEL VARIABLES 

Description Stability Model Variable 


Balance piston damping curvefit coefficients ACB 

Curvefit coefficients for bearing/turbine force curve AK 

Balance piston stiffness curvefit coefficients AKB 

Absolute value of OM AOM 

* AOM to the 0,1, 2, 3, 4 powers AOMP 

Scale factor for Alford cross-coupling coefficient BETA 

Pump bearing damping Cl 

Turbine bearing damping C2 

Modal damping matrix CMAT 

Array of damping matrices at 7 joints CS 

Rotor modal mass integrals GA0.GA1 

Modal mass integral constants 11,12 

Stiffness of backup structure at bearings KBUP,KBUT 

Modal stiffness matrix KMAT 

Array of stiffness matrices at 7 joints KS 

Bearing stiffness KTDB 

Rigid rotor angular velocity OM 

Case mode natural frequencies OMC 

Rotor mode natural frequency OMR 

Mode shapes PHI 

Alford force cross-coupling curve fit coeff, QUALF 

RPM - rad conversion RPRD 


El 





ACBO.ACBl ,ACB2 

AKBO.AKBl ,AKB2 

AKTY.AKTZ 

AKPY(0-4) 

AKPZ(0-4) 

AKTY(0-4) 

AKTZ(0-4) 

BETA 

BKPY.BKPZ 

BKTY.BKTZ 

CPY.CPZ 

CTY.CTZ 

GAMO 

GAM1 

11 

12 

PHICPY.PHICPZ 

PHICTY.PHICTZ 

PHICBX 


ACB[1-3J 
AKB[ 1-3] 

KTDB 
AK[ 1 ; ] 

AK[2;] 

AK[3;] 

AK[4 ; ] 

BETA 

KBUP 

KBUT 

Cl 

C2 

GAO 

GA1 

11 

12 

(M**-.5)PHI[1;2,3;8-17] 
(M**-.5)PHI[2;2,3;8-17] 
(M**-.5)PHI[3;1 ;8-17] 
(M**-.5)PHI[3;2,3;8-17] 

( I 2**-. 5 ) PHI [4; 2, 3 ;8- 17] 
( 12**-. 5) PHI [5; 2, 3 ;8- 17] 


PHICBY.PHIBPZ 

PHICAY.PHICAZ 

PHICS1Y.PHICS1Z 




SSME VARIABLE COMPARISON, con't 
Hybrid Simulation 

PHICS2Y.PHICS2Z 

PHICS3Y,PHICS3Z 

PHIRPY.PHIRPZ 

PHIRTY.PHIRTZ 

PHIRBX 

PHIRBY.PHIRBZ 

PHIRAY.PHIRAZ 

PHIS1Y.PHIS1Z 

PHIS2Y.PHIS2Z 

PHIS3Y.PHIS3Z 

PHIDXA 

QA( 0-4 ) 

SCA(0-2) 

SC 1 ( 0—2) 

SC2( 0-2) 

SC3( 0-2) 

SKA(0-2) 

SKI ( 0-2) 

SK2( 0-2) 

SK3(0-2) 

SQA(0-2) 

SQl(0-2) 

SQ2(0-2) 

SQ3( 0-2 ) 

TWOZETC 

TWOZETR 


PDP Stability Model 

(M**-.5)PHI[6;2,3;8-17] 
(M**-.5)PHI[7;2,3;8-17] 
(M**-.5)PHI[1 ;2,l-7] 
(M**-. 5) PH I [2; 2, 3; 1-7] 
(M**-.5)PHI[3;1 ;l-7] 
(M**-.5)PHI[3;2,3;l-7] 

( 1 2**-. 5) PHI [4; 2, 3; 1-7] 
(1 2**-. 5) PHI [5; 2, 3; 1-7] 
(M**-. 5) PHI [6; 2, 3; 1-7] 
(M**-. 5) PHI [7 ;2 , 3 ; 1-7] 
AOM 

24*QALF( 1-5) 

SC[;4] 

sc[ ; l] 
sec *2] 
sec ; 3] 

SK[;4] 

SK[;1] 

SK[;2] 

SK[ ; 3] 

SQC ;4] 

SQC ; l] 

SQC; 2] 

SQC;3] 

2*ZEC 

2*ZER 


F2 


